Measurement of Neutrino Oscillation by the K2K Experiment 
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We present measurements of disappearance in K2K, the KEK to Kamioka long-baseline neu- 
trino oscillation experiment. One hundred and twelve beam-originated neutrino events are observed 
in the fiducial volume of Super-Kamiokande with an expectation of ISS.ltg'e events without os- 
cillation. A distortion of the energy spectrum is also seen in 58 single-ring muon-like events with 
reconstructed energies. The probability that the observations are explained by the expectation for 
no neutrino oscillation is 0.0015% (4.3(t). In a two flavor oscillation scenario, the allowed Am^ 
region at sin^ 26 = 1 is between 1.9 and 3.5 x 10"'^ eV'^ at the 90 % C.L. with a best-fit value of 
2.8 X 10"^ eV^ 

PACS numbers: 14.60.Pq,13.15.+g,25.30.Pt,95.55.Vj 



I. INTRODUCTION 

The oscillation of neutrinos into other neutrino fla- 
vors is now well established. By using the angle and 
energy distribution of atmospheric neutrinos, the Super- 
Kamiokande collaboration has measured the parameters 
of oscillation and observed the sinusoidal disappearance 
signature predicted by oscillations [l|, Q ■ The K2K col- 
laboration has previously reported evidence of neutrino 
oscillations in a man-made neutrino beam which was di- 
rected 250 km across Japan [l, |3| . 

For neutrinos of a few GeV, the dominant oscillation is 
between and Vr flavor states and two-flavor oscillations 
suffice to describe and analyze the data. In the two- 
flavor neutrino oscillation framework the probability that 
a neutrino of energy E^, with a flavor state will later 
be observed in the v-r flavor eigenstate after traveling a 
distance L in vacuum is: 



P{v^ -> Vr) = sin^ 20 am^ 



1.27ATO2(eV^)L(km) 
^.(GeV) 



, (1) 



where 6 is the mixing angle between the mass eigenstates 
and the flavor eigenstates and Am^ is the difference of 
the squares of the masses of the mass eigenstates. 

The KEK to Kamioka long-baseline neutrino oscilla- 
tion experiment (K2K) [5| uses an accelerator-produced 
beam of nearly pure with a neutrino flight distance of 
250 km to probe the same Am^ region as that explored 
with atmospheric neutrinos. The neutrinos are measured 
first by a suite of detectors located approximately 300 
meters from the proton target and then by the Super- 
Kamiokande (SK) detector 250 km away. The near de- 
tector complex consists of a 1 kiloton water Cherenkov 
detector (IKT) and a fine grained detector system. SK 
is a 50 kiloton water Cherenkov detector, located 1000 m 
underground 



The K2K experiment is designed to measure neutrino 
oscillations using a man-made beam with well controlled 
systematics, complementing and confirming the measure- 
ment made with atmospheric neutrinos. In this paper we 
report a complete description of the observation of neu- 
trino oscillations in the K2K long-baseline experiment, 
and present a measurement of the Am^ and mixing an- 
gle parameters. 

Neutrino oscillation causes both a suppression in the 
total number of events observed at SK and a dis- 
tortion of the measured energy spectrum compared to 
that measured at the production point. Therefore, all 
of the beam-induced neutrino events observed within the 
fiducial volume of SK are used to measure the overall 
suppression of fiux. In addition, in order to study the 
spectral distortion, the subset of these events for which 
the incoming neutrino energy can be reconstructed are 
separately studied. 

If the neutrino interaction which takes place at SK is a 



charged-current (CC) quasi-elastic(QE)(i/^ 



the incoming neutrino energy can be reconstructed using 
two-body kinematics, and the spectral distortion stud- 
ied. At the energy of the K2K experiment typically only 
the muon in this reaction is energetic enough to produce 
Cherenkov light and be detected at SK but kinematics of 
the muon alone are enough to reconstruct the energy for 
these events. 

In order to select the charged-current (CC) quasi- 
elastic (QE) events in the data sample, one-ring events 
identified as a muon (li?/i) are chosen which have a high 
fraction of CC-QE at the K2K energy. For these events, 
the energy of the parent neutrino can be calculated by 
using the observed momentum of the muon, assuming 
QE interactions, and neglecting Fermi momentum: 



mNE^ - 
ruN - Ef, + Pf^ cos 9fj, ' 



(2) 
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FIG. 1: A schematic view of the KEK-PS and neutrino beam 
line and the location of beam line components. The EPl 
neutrino beam line leads protons through a distance of 400 m 
from the EPl-A extraction point to the target station via the 
straight and arc sections. The characters "C" and "S" in the 
figure show the locations of the CT and SPIC installations, 
respectively. The lower-left inset is a magnified view of the 
target station. The production target and a set of horn mag- 
nets are located in the target station. A pion monitor was 
installed on two occasions downstream the horn magnets. 



where mj\f, E^, m^, and 6^ are the nucleon mass, 
muon energy, the muon mass, the muon momentum and 
the scattering angle relative to the neutrino beam direc- 
tion, respectively. 

In this paper, all data taken in K2K between June 
1999 and November 2004 are used to measure the sup- 
pression of events and energy distortion and to measure 
the parameters of oscillation. 



II. NEUTRINO BEAM 
A. K2K neutrino beam and beam monitor 

The accelerator and the neutrino beam line for K2K 
consist of a 12 GeV proton synchrotron (KEK-PS), a 
primary proton transportation line, a hadron production 
target, a set of focusing horn magnets for secondary par- 
ticles, a decay volume, and a beam dump. A schematic 
view of the KEK-PS and neutrino beam line is shown 
in Fig. [TJ In this section, we describe each beam line 
component in order, from upstream to downstream. 



1. Primary proton beam 

Protons are accelerated by the KEK-PS to a kinetic 
energy of 12 GeV. After acceleration, all protons are ex- 
tracted in a single turn to the neutrino beam line. The 
duration of an extraction, or a "spill", is 1.1 fj,sec, which 
contains 9 bunches of protons with a 125 ns time interval 



between them. As shown in Fig.[l] the beam is extracted 
toward the north, bent by 90° toward the direction of 
SK, and transported to the target station. There is a 
final steering magnet just before the target which directs 
the beam to SK at an angle of about 1° downward from 
horizontal. 

The beam intensity is monitored by 13 current trans- 
formers (CTs) installed along the neutrino beam line as 
shown in Fig. [TJ The CTs are used to monitor the beam 
transportation efhciency. The overall transportation ef- 
ficiency along the beam line is about 85%. A CT placed 
just in front of the production target is used to estimate 
the total number of protons delivered to the target. A 
typical beam intensity just before the target is about 
5 X 10^^ protons in a spill. 

In order to measure the profile and the position of the 
beam, 28 segmented plate ionization chambers (SPICs) 
are also installed (Fig. [T]). They are used to steer and 
monitor the beam, while the last two SPICs in front of the 
target are used to estimate the beam size and divergence, 
which is used as an input to our beam Monte Carlo (MC) 
simulation. 



2. Hadron production target and horn magnets 

A hadron production target and a set of horn magnets 
are placed in the target station. Protons hit the target 
and a number of secondary particles are generated at the 
production target. Two toroidal magnetic horns are em- 
ployed to focus positively charged particles, mainly 7r"'"'s, 
in the forward direction by the magnetic field. A typi- 
cal focusing of transverse momentum by the horn mag- 
nets is about 100 MeV/c per meter. The momenta of 
focused pions are around 2—3 GeV/c, which corresponds 
to about 1.0—1.5 GeV of energy for those neutrinos de- 
caying in the forward direction. According to our Monte 
Carlo simulation, the flux of neutrinos above 0.5 GeV is 
22 times greater with horn magnets with 250 kA current 
than without the horn current. 

A schematic view of the horn magnets is shown in 
Fig. [21 The dimensions of the first horn are 0.70 m in 
diameter and 2.37 m in length, while those of the second 
horn are 1.65 m in diameter and 2.76 m in length. Both 
horns are cylindrically symmetric in shape. The produc- 
tion target, a rod of a length of 66 cm and diameter of 
3 cm, made of aluminum alloy 606 1-T, is embedded in- 
side the first horn. The target diameter was 2 cm in June 
1999 and was changed to 3 cm in November 1999 for im- 
proved mechanical strength. The target also plays the 
role of inner conductor of the first horn, making a strong 
magnetic field inside the horn to achieve high focusing ef- 
ficiency. The second horn is located 10.5 m downstream 
from the first horn, playing the role of a refiector, which 
re-focuses over-bent low energy pions, and in addition 
further focuses under-bent high energy pions. 

Pulsed current with a duration of 2 msec and an ampli- 
tude of 250 kA (200 kA in June 1999) is supphed by four 
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2nd Magnetic Horn 



1st Ivlagnetic Horn 



Production Target 




FIG. 2: Schematic view of the two horn magnets. An elec- 
trical current of 250 kA is supplied to both horns, creating 
a toroidal magnetic field inside the horns. The production 
target, an aluminum rod of 66 cm in length and 3 cm in di- 
ameter, is embedded inside the first horn magnet, which also 
plays the role of inner conductor of the horn. The second 
horn is located 10.5 m downstream of the first horn. 



current feeders to each horn. The peaking time of the 
current is adjusted to match the beam timing. The max- 
imum magnetic field in the horn is 33 kG at the surface 
of the target rod with 3 cm diameter target and 250 kA 
horn current. 

The values of the current supplied to the horn mag- 
nets are read out by CTs put in between current feed- 
ers and recorded by a flash analog-to-digital converter 
(FADC) on a spill-by-spill basis. Overall current and 
current balance between feeders are monitored to select 
good beam spills. The magnetic field inside the proto- 
type of the first horn was measured using pickup coils; 
results showed that the radial distribution of the field 
was in agreement with the design distribution and the 
azimuthal symmetry was confirmed to within a measure- 
ment error of 15%. Detailed descriptions of the horn 
magnets are found in 7, 8, !|. 

A pion monitor (PIMON) was installed on two occa- 
sions just downstream of the horn magnets, as shown in 
Fig. [U in order to measure the momentum and angular 
distributions of pions coming through the horn magnets. 
The PIMON wiU be described in detail later in Sec. IVl 



3. Decay volume, beam dump, and muon monitors 

The positive pions focused by the horn magnets go into 
a 200 m long decay volume which starts 19 m downstream 
of the production target, where the tt"*" decay: 7r+ 
fj,'^ . The decay volume is cylindrical in shape and is 
separated into three sections with different dimensions. 
The diameters of the pipe are 1.5 m, 2 m, and 3 m in 
the first 10 m, the following 90 m, and the remaining 
100 m sections, respectively. The decay volume is filled 
with helium gas of 1 atm (rather than air) to reduce the 
loss of pions by absorption and to avoid uncontrollable 
pion production in the gas. The beam dump is located 
at the end of the decay volume to absorb all the particles 
except for neutrinos. It consists of 3.5 m thick iron, 2 m 
thick concrete, and a region of soil about 60 m long. 

There is a pit called the "muon-pit" just downstream 
of the iron and concrete shields. Muons with momentum 
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FIG. 3: A schematic view inside the muon-pit. An ionization 
chamber (ICH) and an array of silicon pad detectors (SPD) 
are located inside the muon-pit. 



greater than 5.5 GeV/c can reach the muon-pit. The flux 
at the pit is roughly 10'* muons/cm^ /spill. The parent 
particles of both muons and neutrinos are pions, so the 
profile center of muons corresponds to that of neutrinos. 
A change in the beam direction by 3 mrad corresponds 
to a change in the neutrino fiux and spectrum at SK of 
about 1%, and hence it must be controlled and monitored 
to be within 3 mrad. Fig. [3] shows a schematic view 
inside the pit. Two detectors (MUMONs) are installed 
in it: one is an ionization chamber (ICH) and the other 
is an array of silicon pad detectors (SPD). The purpose 
of these detectors is to measure the profile and intensity 
of muons penetrating the shields on spill-by-spill basis. 

An ICH is a segmented plate chamber with a size of 
190 cm (horizontal) x 175 cm (vertical). It consists of six 
modules of size 60 cm x 90 cm, 3 modules in the horizon- 
tal direction and 2 modules in the vertical direction. The 
gap between modules is 25 cm in horizontal and 15 cm 
in vertical (Fig. [3]). The corresponding strip lines of ad- 
joining modules are electrically connected over the gaps 
to make long strip lines of length of ^-^180 cm. There are 
36 horizontal readout channels and 32 vertical channels. 
The channel-to-channel uniformity is calibrated by mov- 
ing ICH horizontally and vertically [lOi] assuming stabil- 
ity of the muon beam. The relative gain of the channels 
has been stable within an accuracy of several percent. 

Two types of SPDs are used: one is a small SPD which 
has a sensitive area of 1 cm x 2 cm with a depletion layer 
thickness of 300 /im, and the other is a large SPD which 
has a sensitive area of 3.4 cm x 3.05 cm with a deple- 
tion layer thickness of 375 fj,m. Seventeen small SPDs 
are arranged along the horizontal and the vertical axes 
at 35 cm intervals while nine large SPDs are in diago- 
nal arrays at 74.2 cm intervals. The sensitivity of each 
small SPD was measured using an LED light source at 
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a test bench and it was found that all the small SPDs 
agree within 6% The sensitivity difference between 
the large SPDs was measured using the muon beam at 
the muon-pit. All the large SPDs were aligned along the 
beam axis simultaneously and the output charge from 
each SPD was compared to obtain the relative gain fac- 
tor. The gain factors have an uncertainty of 10% due to 
the z-dependence of the muon beam intensity [l0| . 



B. Summary of beam operation 

The construction of neutrino beam line was completed 
early in 1999 and beam commissioning started in March 
1999. The beam line and all the components were con- 
structed and aligned within an accuracy of 0.1 mrad with 
respect to a nominal beam axis which was determined 
based on the results of a global positioning system (GPS) 
surve y ac curate to 0.01 mrad between KEK and Kamioka 
sites [Il[. In June 1999, the neutrino beam and detec- 
tors were ready to start data-taking for physics. We took 
data on and off over the period from June 1999 to Novem- 
ber 2004, which is divided into five subperiods according 
to different experimental configurations: June 1999 (la), 
November 1999 to July 2001 (lb), December 2002 to June 
2003 (Ila), October 2003 to February 2004 (lib), and Oc- 
tober 2004 to November 2004 (lie) . The horn current was 
200 kA (250 kA) and the diameter of the production tar- 
get was 2 cm (3 cm) in the la (other) period. The SK 
PMTs were full density for la and lb, but were half den- 
sity for Ila, lib and lie. There was a lead-glass calorime- 
ter (LG) installed in between a scintillating- fiber/water- 
target tracker (SciFi) and a muon range detector (MRD) 
during the la and lb periods; it was replaced by a to- 
tally active fine-segmented scintillator tracker (SciBar) 
for Ila, lib and lie. Only the first four layers of the 
SciBar detector were installed for Ila while it was in its 
full configuration for lib and lie. Furthermore, the wa- 
ter target in the SciFi was replaced by aluminum rods 
during lie. The different experimental configurations for 
the different periods are briefly summarized in Table |T1 

The number of protons delivered to the target is sum- 
marized in Table [H and shown as a function of time 
in Fig. 21 Among the delivered spills, spills which sat- 
isfy the following criteria are used for the physics anal- 
ysis: (1) beam spills with normal machine status. Spills 
during machine studies, beam tuning, and several beam 
studies are discarded. (2) Beam spills with no trouble 
in the beam components and data acquisition systems. 
(3) Beam spills with the proton intensity greater than 
1 X 10^^ protons. (4) Beam spills with the horn cur- 
rent greater than 240 kA (190 kA) for the period other 
than la (for the la period). The number of protons on 
target (POT) for the physics analysis is summarized in 
Table U as well as the total number of protons delivered. 
In total, 1.049 x 10^" protons were delivered to the pro- 
duction target while 0.922 x 10^" POT are used in our 
physics analysis. 
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FIG. 4: The number of protons delivered to the production 
target in the period from March 1999 to November 2004. The 
horizontal axis corresponds to the date. The upper figure 
shows the total number of protons on target (POT) accumu- 
lated since March 1999, and the lower figure shows the POT 
per spill averaged in a day. In total, 104.90 x lO'^^ protons 
were delivered during the entire period including beam com- 
missioning and tuning periods. 



During these periods, the direction of the neutrino 
beam was monitored by MUMON in the muon-pit. Fig. [5] 
shows the stability of the center of the muon profile mea- 
sured by the ionization chamber (ICH) in MUMON. The 
beam was pointed to the direction of SK within ±1 mrad 
during the whole run period, so that the neutrino flux 
and spectrum at SK was stable within much better than 
1%. 



C. K2K neutrino beam simulation 

We use a neutrino beam Monte Carlo (beam MC) sim- 
ulation program to study our neutrino beam properties. 
The beam line geometry is implemented in GEANT [l3| 
and particles are tracked in materials until they decay 
into neutrinos or are absorbed in the material. The tracks 
of neutrinos are extrapolated along a straight line to the 
near detector (ND) and Super-Kamiokande (SK) and the 
fluxes and the energy spectrum at these locations are de- 
termined. 

In the simulation program, protons with a kinetic en- 
ergy of 12 GeV are injected into the aluminum produc- 
tion target. The proflle and divergence are assumed to 
be Gaussian-like and the values measured by two SPICs 
in front of the target are used as inputs. An empirical 
formula for the differential cross-section by J. R. Sanford 
and C. L. Wang [13, [3] is used to simulate the primary 
hadron production in the target. The Sanford- Wang for- 
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TABLE I: Summary of the number of protons on target and the experimental configuration for each running period. The row 
labeled "LG/SciBar configuration" indicates the detector installed between the SciFi and MRD detectors. For the row "SK 
configuration", "SK-I" refers to the configuration with full PMT density while "SK-II" refers to that with half density. See 
the text for more detailed description of the experimental configurations. The delivered POT shown in the table includes the 
beam delivered during commissioning and beam tuning work before the physics runs. 



Periods 


la 


lb 


Ila 


lib 


lie 


total 


Jun.'99 


Nov.'99-JuI.'Ol 


Dec.'02-Jun.'03 


Oct.'03-Feb.'04 


Oct.'04-Nov.'04 


Delivered POT (xlO'**) 


6.21 


49.85 


24.91 


20.15 


3.78 


104.90 


POT for analysis (xlO^*) 


3.10 


44.83 


22.57 


18.61 


3.12 


92.23 


Horn current 


200 kA 


250 kA 


250 kA 


250 kA 
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FIG. 5: Stability of the center of muon profile measured by 
the ionization chamber (ICH) in MUMON. The upper figure 
shows the profile center of the horizontal direction and the 
lower figure shows that of the vertical direction. In each fig- 
ure, the beam direction to SK measured by GPS and ±1 mrad 
off the center are indicated by horizontal lines. The data 
shown here are after good beam selection. 



mula is expressed as following: 
dildp \ pb J 



(3) 



exp 



Pi 



Cq9{p- Cj Pb cos 



where (Pa/dfldp is the double differential cross section of 
particle production per interacting proton in the unit of 
mb sr^^ (GcV/c)~^, 9 is the angle between the secondary 
particle and the beam axis in the laboratory frame, p and 
Pb are the momenta of the secondary particle and the 



TABLE II: The fitted parameters, d's, in the Sanford-Wang 
formula for the production of positively charged pious in the 
Cho-CERN compilation and for the HARP results [Hi. The 
target nucleus is beryllium in Cho-CERN compilation while 
it is aluminum in the HARP results. The values in the table 
are before the nuclear scaling is applied. 





Ci 


C2 


Ca 


Ci 


Cs 


Ce 


Ct 


Cs 


HARP 
Cho-CERN 


440 
238 


0.85 
1.01 


5.1 
2.26 


1.78 
2.45 


1.78 
2.12 


4.43 
5.66 


0.14 
0.14 


35.7 
27.3 



incident proton, respectively. The C^'s are parameters 
fitted to existing hadron production data. For the pro- 
duction of positively charged pions, we use as a reference 
model the Q's obtained from a fit designated the "Cho- 
CERN compilation" , in which the data used in the com- 
pilation mainly come from the measurement of proton- 
beryllium interactions performed by Cho et al. [lH . The 
values for C^'s are shown in Table [Til A nuclear rescal- 
ing is then applied to convert the pion production cross 
section on beryllium to that on aluminum. The scaling 
factor, w, is defined as 



(4) 



where Aai and are atomic masses for aluminum and 
beryllium, respectively, and an index a^xp) is expressed 
as 



aixp) = 0.74 + xf{-0.55 + 0.26xf) 



(5) 



as a function of the Feynman x variable, xp- 

Negatively charged pions and charged and neutral 
kaons are generated as well as positively charged pions 
using the same Sanford-Wang formula with different sets 
of Ci's. For negative pion production, t he p arameters in 
are used, while those described in [l^ are used for 
the kaon production. 

Generated secondary particles are tracked by GEANT 
with the GCALOR/FLUKA [H, [H, hadron model 
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FIG. 6: The energy spectrum for each type of neutrino at 
ND (left) and SK (right) estimated by the beam MC simu- 
lation. The neutrino beam is 97.3% (97.9%) pure muon neu- 
trino with contaminations of v^jvp, ~ 0.013 (0.009), V^/Uf_i ~ 
0.015 (0.012), and Ve/vt. ~ 1.8 x 10"" (2.2 x 10"*) at ND 
(SK). 



SciFi Detector 




FIG. 7: The schematic view of the near neutrino detectors 
for K2K-IIb period. In K2K-I, the Lead-Glass calorimeter 
was located at the position of the SciBar detector. 



through the two horn magnets and the decay volume until 
they decay into neutrinos or are absorbed in materials. 

Since GEANT treats different types of neutrinos iden- 
tically, we use a custom-made simulation program to 
treat properly the type of neutrinos emitted by particle 
decays. Charged pions are treated so that they decay into 
muon and neutrino (7r+ — > fj.'^v^^TT^ — s- fi~T'^^ called 
TT^j) with branching fraction of 100%. The kaon de- 
cays considered in our simulation are so- 

called ii:;^2> ^^3° 

and K^^^* decays. Their branching ratios are taken from 
the Particle Data Group [2l[. Other decays are ignored. 
Neutrinos from Kg are ignored since the branching ratio 
for Kg decaying to neutrinos is quite small. The Dalitz 
plot density of V—A theory [U |23| is employed prop- 
erly in Ke3 decays. Muons are considered to decay via 
//^ Vei'^'e) called with 100% branch- 

ing fraction. The energy and angular distributions of the 
muon antineutrino (neutrino) and the electron neutrino 
(antineutrino) emitted from a positive (negative) muon 
are calculated according to Michel spectra of V—A the- 
ory (2^ , where the polarization of the muon is taken into 
account. 

The produced neutrinos are extrapolated to the ND 
and SK according to a straight line and the energy and 
position of the neutrinos entering the ND and SK are 
recorded and used in our later simulations for neutrino 
interaction and detector simulators. 

The composition of the neutrino beam is dominated 
by muon neutrinos since the horn magnets mainly fo- 
cus the positive pions. Figure [5] shows the energy spec- 
tra of each type of neutrino at ND and SK estimated 
by the beam MC simulation. About 97.3% (97.9%) 
of neutrinos at ND (SK) are muon neutrinos decayed 
from positive pions, and the beam is contaminated with 
a small fraction of neutrinos other than muon neutri- 
nos; Ve/v^, - 0.013 (0.009), Vf,/Vf, ~ 0.015 (0.012), and 
V^/v^ - 1.8 X 10-4 (2.2 X 10-"*) at ND (SK). The va- 
lidity of our beam MC simulation has been confirmed by 



both the HARP experiment and PIMON measurements, 
which will be described in detail in Sec. |Vl 



III. NEUTRINO DETECTORS 

A near neutrino detector system (ND) is located 300 m 
downstream from the proton target. The primary pur- 
pose of the ND is to measure the direction, flux, and the 
energy spectrum of neutrinos at KEK before they oscil- 
late. The schematic view of the ND during the K2K- 
Ilb period is shown in Fig. [71 The ND is comprised 
of two detector systems; a one kiloton water Cherenkov 
detector (IKT) and a fine-grained detector (FGD) sys- 
tem. The FGD consists of a scintillating-fiber/water- 
target tracker (SciFi), a Lead-Glass calorimeter (LG) 
in K2K-I period, a totally active fine-segmented scin- 
tillator tracker (SciBar) in K2K-IIb and K2K-IIc peri- 
ods, and a muon range detector (MRD). The far detec- 
tor is the 50 kiloton water Cherenkov detector, Super- 
Kamiokande (SK), which is located 250 km away from 
KEK and 1000 m (2700 m water equivalent) below the 
peak of Mt. Ikeno-yama in Gifu prefecture. 



A. 1 kiloton water Cherenkov detector 

A one kiloton water Cherenkov detector (IKT) is lo- 
cated in the experimental hall at KEK as the upstream 
detector. The IKT detector is a miniature version of 
SK, and uses the same neutrino interaction target mate- 
rial and instrumentation. The primary role of the IKT 
detector is to measure the j/^ interaction rate and the 
energy spectrum. The IKT detector also provides a high 
statistics measurement of neutrino- water interactions. 

The cylindrical tank, 10.8 m in diameter and 10.8 m 
in height, holds approximately 1000 tons of pure water. 
The center of the water tank is 294 m downstream of the 
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pion production target. The water tank is optically sepa- 
rated into the inner detector (ID) and the outer detector 
(OD) by opaque black sheets and reflective Tyvek® (a 
material manufactured by DuPont) sheets. The ID of the 
IKT detector is a cylinder of 8.6 m in diameter and 8.6 m 
in height. This volume is viewed by 680 photomultiplier 
tubes (PMTs) of 50 cm diameter facing inward to detect 
Cherenkov light from neutrino events. The PMTs and 
their arrangement are identical to those of SK; 70 cm 
spacing between PMTs gives a 40% photocathode cov- 
erage. The fiducial volume used for selecting neutrino 
events in the IKT is defined as a 25 ton cylindrical re- 
gion with a diameter of 4 m and a length of 2 m oriented 
along the beam axis. The OD covers the upstream third 
of the barrel wall and the whole of the bottom wall. The 
OD volume is viewed by 68 PMTs of 20 cm diameter, 
facing outward to veto the incoming particles. The OD 
is also used to trigger through-going/stopping cosmic ray 
muon events for detector calibrations. 

To compensate for the geomagnetic field which affects 
the PMT response, nine horizontal Helmholtz coils and 
seven vertical Helmholtz coils are arranged surrounding 
the water tank. The water purification system for the 
IKT detector circulates about 20 tons/hour of water. 
The electrical resistance (~10 Mfi/cm) and water tem- 
perature (^11 °C) are kept constant by the system. 

The IKT detector data acquisition (DAQ) system is 
similar to that of SK. The signal from each PMT is pro- 
cessed using custom electronics modules called ATMs, 
which were developed for the SK experiment and are 
used to record digitized charge and timing information 
for each PMT hit over a threshold of about 1/4 photo- 
electrons. The DAQ trigger threshold is about 40 PMT 
hits within a 200 nsec time window in a 1.2 /isec beam 
spill gate, where the beam spill gate is issued to all near 
detectors, synchronized with the beam timing, by the ac- 
celerator. The 40 hit threshold is roughly equivalent to 
the signal of a 6 MeV electron. The pulse shape of the 
analog sum of aU 680 PMTs' signals (PMTSUM) is also 
recorded for every beam spill by a 500 MHz flash analog 
to digital converter (FADC) which enables us to iden- 
tify multiple interactions in a spill gate. We determine 
the number of interactions in each spill by counting the 
peaks in PMTSUM greater than a threshold equivalent 
to a 100 MeV electron signal. 

The physical parameters of an event in the IKT detec- 
tor such as the vertex position, the number of Cherenkov 
rings, particle types and momenta are determined us- 
ing the same algorithms as in SK [T]. First, the vertex 
position of an event is determined from the PMT tim- 
ing information. With knowledge of the vertex position, 
the number of Cherenkov rings and their directions are 
determined by a maximum-likelihood procedure. Each 
ring is then classified as e-like, representing a showering 
particle (e^ , 7) , or ^-like, representing a non-showering 
particle (/x^,7r*), using its ring pattern and Cherenkov 
opening angle. On the basis of this particle type infor- 
mation, the vertex position of a single-ring event is fur- 



ther refined. The momentum corresponding to each ring 
is determined from the Cherenkov light intensity. Fully 
contained (FC) neutrino events, which deposit all of their 
Cherenkov light inside the inner detector, are selected by 
requiring the maximum number of photoelectrons on a 
single PMT at the exit direction of the most energetic 
particle to be less than 200. The events with the maxi- 
mum number of photoelectrons greater than 200 are iden- 
tified as a partially contained (PC) event. This criterion 
is used because a muon passing through the wall produces 
a lot of light in the nearest PMTs. 

The reconstruction quality, especially the vertex posi- 
tion and angular resolution, are estimated with a MC 
simulation. The vertex resolution is estimated to be 
14.7 cm for FC single-ring events and 12.5 cm for PC 
single-ring events, while those for multi-ring FC and PC 
events are 39.2 cm and 34.2 cm, respectively. The angu- 
lar resolution for single-ring CC-QE events is estimated 
to be 1.05° for FC events and 0.84° for PC events. As 
for the capability of the particle identification, 0.3% of 
muon neutrino CC quasi-elastic events with a single ring 
are misidentified as e-like while 3.3% of electron neutrino 
CC quasi-elastic events with a single ring are misidenti- 
fied as /i-like. The momentum resolution for muons is 
estimated to be 2.0-2.5% in the whole momentum range 
of the IKT. 

The gain and timing of each PMT are calibrated us- 
ing a Xe lamp and a N2 laser as light sources, respec- 
tively. The absorption and scattering coefficients of water 
are measured using laser calibration, and the coefficients 
in the detector simulation are further tuned to repro- 
duce the observed charge patterns of cosmic ray muon 
events. The energy scale is calibrated and checked by 
cosmic ray muons with their decay electrons and neutral 
current 7r°s produced by the K2K neutrino beam. The 
absolute energy scale uncertainty is ^4% while the verti- 
cal/horizontal detector asymmetry of the energy scale is 
1.7%. The energy scale is stable within about 1% from 
2000 to 2004. 

The performance of vertex reconstruction is experi- 
mentally studied by special cosmic ray muon data uti- 
lizing a PVC pipe with scintillating strips at each end 
inserted vertically into the tank. Cosmic ray muons go- 
ing through the pipe emulate the neutrino-induced muons 
whose vertex position is defined at the bottom end of the 
pipe. This study demonstrates that the vertex recon- 
struction works as well as we expected from the Monte 
Carlo simulation. We find a vertex bias difference be- 
tween data and MC simulation of less than 4 cm for both 
FC and PC events. 



B. Scintillating Fiber detector 

The scintillating fiber (SciFi) detector is a 6 ton track- 
ing detector with integral water target layers. Details of 
the design and performance of the detector are described 
in Refs.|23l [23 |. The SciFi detector is used to measure 
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the neutrino spectrum, and to reconstruct with high res- 
olution the charged particle tracks produced in neutrino 
interactions. It can estimate the rates for quasi-elastic 
and inelastic interactions and is sensitive to higher en- 
ergy events, and hence has complementary capabilities 
to the IKT detector. The SciFi detector has been in 
stable operation since 1999 when the first K2K neutrino 
beam was delivered. 

The SciFi detector consists of 20 layers of 2. 6m x 2.6m 
tracking modules, placed 9 cm apart. Each layer contains 
a double layer of sheets of scintillating fibers arranged, 
one each, in the horizontal and vertical directions; each 
sheet is itself two fibers thick. The diameter of each fiber 
is 0.692 mm. In between the fiber modules, there are 19 
layers of water target contained in extruded aluminum 
tanks. The water level was monitored; it has stayed con- 
stant within 1% throughout the experiment, except for a 
few tanks which drained following an earthquake. This 
monitoring, as well as measurements when the tanks were 
filled and later drained, give a fiducial mass of 5590 kg 
with 1% accuracy. The fiducial mass fractions are 0.700 
H2O, 0.218 Al, and 0.082 HC (± 0.004). 

The fiber sheets are coupled to an image intensifier 
tube (IIT) with a CCD readout system. The relative 
position between the fibers and the CCD coordinate sys- 
tem is monitored periodically by illuminating every 10th 
or 20th fiber with an electro-luminescent plate placed at 
the edge of each fiber sheet. In addition, cosmic-rays 
were used to monitor the gain of the system on a weekly 
basis. 

Hit fibers are extracted using the CCD images. The 
raw data consists of hit CCD pixels and their digitized 
brightness. Neighboring hit pixels are grouped to make 
a pixel cluster. Those clusters are then combined and 
matched to the location of specific scintillating fibers. 
The efficiency to identify a fiber through which a charged 
particle passed is estimated using cosmic ray muons to 
be about 95%, but closer to 90% at angles within 30 
degrees of the beam. After hit fibers are reconstructed, 
tracks with three or more hit layers are reconstructed 
using conventional fitting techniques. The efficiency to 
find a track is also estimated using cosmic ray muons, and 
is ~ 70% for tracks with length of three layers, ~ 87% 
for four layers, and approaches 100% for longer tracks. 

Surrounding the SciFi are two plastic scintillator ho- 
doscope systems. One is placed downstream of SciFi and 
gives track timing and position information. It also serves 
as a pre-shower detector for the Lead Glass calorimeter. 
The other is upstream of SciFi and is used to veto muons 
and other particles from the beam, primarily from neu- 
trino interactions in the upstream IKT detector, but also 
from cosmic rays. 

The downstream system consists of 40 scintillator units 
placed one upon another having a total height of 4 m. 
Each unit is made of a plastic scintillator 466 cm long, 
10.4 cm high, and 4 cm thick. A PMT is attached to each 
end of the scintillator. The horizontal position of the 
charged particle can be calculated with 5 cm resolution 



from the timing information read out by the both end 
PMT's. The upstream veto wall is similar, but pairs of 
scintillators are joined together by optical cement and 
share a single light guide for each PMT. Thus there are 
fewer readout channels and the vertical resolution is twice 
as coarse, but the hodoscope covers the same total area as 
the one downstream. The charge and timing information 
from each of the 120 total PMT's are recorded. The 
energy deposit measured in the downstream hodoscope is 
used to select electron neutrino events as described later. 
The energy resolution of these hodoscopes is estimated 
using cosmic ray muons to be 7.4% for minimum ionizing 
particles. 

A more detailed description of the hodoscope system 
can be found in [25] . 



C. Scintillating Bar detector 

The SciBar detector [2^ was constructed as an upgrade 
of the near detector system. The purposes of the SciBar 
detector are to measure the neutrino energy spectrum 
and to study the neutrino interaction with high detec- 
tion efficiency for low momentum particles. The main 
part of the SciBar detector consists of an array of plastic 
scintillator strips. Its totally active and finely segmented 
design allows us to detect all the charged particles pro- 
duced in a neutrino interaction. 

We use extruded scintillator strips produced by 
FNAL j23|. The dimensions of a strip are 1.3 cm thick, 
2.5 cm wide, and 300 cm long. In total, 14,848 scintilla- 
tor strips are arranged in 64 layers of alternating vertical 
and horizontal planes. The dimension of the detector is 
3 X 3 X 1.7 m'^ providing the total weight of about 15 tons. 

The scintillation light is guided to multi-anode PMTs 
by wavelength shifting fibers inserted into the holes of 
scintillator strips. Sixty-four wavelength shifting fibers 
are bundled together and glued to an attachment to be 
precisely coupled between fibers and the photo cathode 
of the multi-anode PMT. Both charge and timing of the 
PMT outputs are recorded using custom-made electron- 
ics [l^. The noise level and the timing resolution for 
minimum-ionizing particle signal are about 0.3 photo- 
electrons and 1.3 nsec, respectively. 

The gain of all multi-anode PMT channels was mea- 
sured at a test bench prior to the installation. In or- 
der to monitor and correct gain drift during operation, 
the SciBar is equipped with a gain calibration system 
using LED [29]. The gain stability is monitored with 
precision better than 1%. Cosmic-ray data are collected 
between beam spills to calibrate the multi-anode PMT 
gain and scintillator light yield in-situ. The light yield 
has been stable within 1 % during operation. The light 
attenuation length of the wavelength shifting fiber is also 
measured with cosmic ray muons. It is confirmed to be 
consistent with the test bench measurement done prior 
to the installation. 

An electromagnetic calorimeter (EC) is installed down- 
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stream the tracker part of SciBar to study the amount 
of the electron neutrino contamination in the beam and 
tt'^ production in neutrino interactions. The calorimeter 
is made of bars of dimensions 262 x 8 x 4 cm'^. The bars, 
a sandwich of lead and scintillating fibers, were origi- 
nally built for the "spaghetti" calorimeter of the CHO- 
RUS neutrino experiment at CERN [30*1 . Each bar is read 
out by two PMTs per side. In the SciBar-EC, 32 bars 
are assembled to form a plane of vertical elements, fol- 
lowed by a plane of 30 horizontal bars. The two planes, 
each 4 cm thick, cover an area of 270 x 262 cm^ and 
262 X 250 cm^, respectively. The EC adds eleven ra- 
diation lengths to the tracker part which has about four 
radiation lengths. The response linearity of the EC is un- 
derstood to b e better t han 10%. The energy resolution is 
about 14%/-y/i5[GeV] as measured with a test beam 

To reconstruct neutrino events, hit scintillator strips 
in SciBar with more than or equal to two photoelectrons 
(corresponding to about 0.2 MeV) are selected. Charged 
particles are reconstructed by looking for track projec- 
tions in each of two dimensional view (x-z and y-z). using 
a cellular automaton algorithm [3l|. Then, track candi- 
dates in two views are combined based on matching of the 
track edges in z direction and timing information. Recon- 
structed tracks are required to have hits in more than or 
equal to three consecutive layers. The minimum length 
of reconstructible track is, therefore, 8 cm, which is corre- 
sponding to 450 MeV/c for protons. The reconstruction 
efficiency for an isolated track longer than 10 cm is 99%. 



D. Muon range detector 

The muon range detector (MRD) [s^l has two pur- 
poses. One is to monitor the stability of the neutrino 
beam direction, profile and spectrum by measuring the 
energy, angle and production point of muons produced 
by charged-current neutrino interaction by utilizing its 
huge mass of the iron as the target. The other is to iden- 
tify the muons produced in the upstream detectors and 
to measure their energy and angle with combination of 
other fine grain detectors. This enables us to measure 
the energy of the incident neutrino. 

MRD consists of 12 layers of iron absorber sandwiched 
in between 13 sets of vertical and horizontal drift-tube 
layers. The size of a layer is approximately 7.6 mx 7.6 m. 
In order to have a good energy resolution for the whole 
energy region, the upstream four iron plates are 10 cm 
thick while the downstream eight plates are 20 cm thick. 
The total iron thickness is 2.00 m covering the muon 
energy up to 2.8 GeV. MRD has 6,632 drift tubes, each 
of which is made of aluminum with a cross section of 
5 cmx 7 cm. PIO gas (Ar : CH4 = 90% : 10%) is suppfied 
to all the tubes. The maximum drift time in a tube is 
about 1 /isec. The drift time is digitized by 20 MHz 6-bit 
TDCs. The total weight of iron is 864 tons and the total 
mass of MRD including the aluminum drift tubes is 915 
tons. 



A conventional track finding algorithm is employed to 
reconstruct tracks from hits. The track finding efficiency 
is 66%, 95% and 97.5% for tracks with one, two and 
three traversed iron plate(s), respectively, and it goes up 
to 99% for longer tracks. The range of track is estimated 
using the path length of the reconstructed track in iron. 

Accurate knowledge of the iron-plate weight is neces- 
sary for the measurements of both neutrino interaction 
rate and track range. Relative thickness of each plate was 
studied by comparing the event rate using the neutrino 
beam data. Also, the density was measured directly us- 
ing a sample of the same iron. Combining these studies, 
we quote the weight of the iron plates with an accuracy 
of 1%. The relation between the muon energy and the 
muon range in iron was calculated using a GEANT based 
Monte Carlo code. There is at maximum 1.7% difference 
in the muon range among various calculations. We quote 
the error on energy scale in the range measurement to be 
2.7% by linearly adding these two errors. 

The energy acceptance and resolutions of the MRD 
were studied by a Monte Carlo simulation. The accep- 
tance is ranging from 0.3 GeV to 2.8 GeV while the res- 
olution is 0.12 GeV for forward-going muons. The track 
angular resolution is about 5 degrees and the resolution 
of the vertex point perpendicular to the beam direction 
is about 2 cm. 



E. Lead glass calorimeter 

The Lead Glass (LG) calorimeter was located between 
SciFi and MRD in K2K-I period. The purpose of LG 
is to distinguish electrons from muons by measuring the 
energy deposit. The LG calorimeter is made up of 600 
cells. A LG cell of approximately 12 cmx 12 cmx 34 cm is 
viewed by 3 inch-in-diameter PMT(Hamamatsu, R1652) 
through a light guide cylinder made also by lead glass. 
This LG calorimeter was once used in the TOPAZ exper- 
iment ^33i] and reused for the K2K experiment. 

The LG detector system reads out only the charge in- 
formation for each cell. The absolute energy scale of 9 
standard LG cells out of 600 were calibrated prior to 
installation by using an electron beam from the elec- 
tron synchrotron with the energy range from 50 MeV 
to 1.1 GeV. The resolution was estimated by this pre- 
calibration to be 10% at 1 GeV. Position dependence for 
the energy resolution were also measured to be 4%. The 
other LG cells were relatively calibrated to the standard 
cells by cosmic-ray muons. 

Responses for muons were also calibrated by using 
cosmic-ray muons at KEK prior to installation. The rel- 
ative peak pulse height for PMTs was adjusted to each 
other within 2%. The responses for charged pions were 
checked at different momenta (0.3—2.0 GeV/c) by using 
the KEK test beam, confirmed to be in good agreement 
with the expectation by an MC simulation. 
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F. Super-Kamiokande 

The far detector of the K2K experiment is Super- 
Kamiokande, which is located in the Kamioka Observa- 
tory, operated by the Institute for Cosmic Ray Research, 
University of Tokyo. The SK detector is a cylindrically 
shaped water Cherenkov detector which is 41 m in height, 
39 m in diameter and has a total mass of 50 kilotons 
of water. The water tank is optically separated into a 
cylindrically-shaped inner detector (ID) and outer de- 
tector (OD) by opaque black sheets and Tyvek® sheets 
attached to a supporting structure. The ID is viewed by 
11,146 20-inch PMTs facing inward covering 40% of the 
ID surface from June 1999 to 2001 (called SK-I and K2K- 
I), while it is viewed by 5,182 PMTs enclosed in a fiber 
reinforced plastic and sealed with acrylic covers on their 
front surface, covering 19% of the ID surface from De- 
cember 2002 (SK-II and K2K-II). The transparency and 
the reflection of these covers in water are 97% and 1%, re- 
spectively. In the OD region, outward-facing 1,885 8- inch 
PMTs are attached to the outer side of the supporting 
structure. The performance of OD PMTs is improved 
in SK-II. The fiducial volume is defined to be a cylinder 
whose surface is 2 m away from the ID wall providing a 
fiducial mass of 22.5 kilotons. Details of the detector per- 
formance and systematic uncertainties in SK-I are writ- 
ten in P, Q . For SK-II, these quantities are estimated 
using similar methods as used in SK-I. Momentum res- 
olution for SK-II is slightly worse than SK-I; 2.4% and 
3.6% for 1 GeV/c muons in SK-I and SK-II, respectively. 
This is because the number of ID PMTs in SK-II is about 
a half of SK-I. However, the performance of the vertex 
reconstruction, the ring counting, and the particle iden- 
tification in SK-II are almost the same as in SK-I. The 
purity of the QE interaction in 1-ring ^-like events is 
58%. The uncertainty in the energy scale is estimated to 
be 2.0% for SK-I and 2.1% for SK-II. 

In this long baseline experiment, timing information 
is used to distinguish between beam neutrino events and 
cosmic ray induced background events in the SK detec- 
tor. The GPS is used to synchronize the timing of the 
beam spill between KEK and SK. At both sites are a free 
running 50 MHz (32-bit) local time counter connected to 
a GPS receiver and an event trigger (at Super-K) or the 
beam spill trigger (at KEK). At first, a quartz oscillator 
was used with good results, and later oscillator drift was 
improved further with a rubidium clock. This counter is 
synchronized using the one pulse-per-second signal from 
the GPS. In this way, events can be synchronized within 
approximately 50 ns, after compensating for oscillator 
drift. This is confirmed by comparing a second, inde- 
pendent timing system at each site which gives the same 
result as the primary system within 35 ns 99% of the 
time. As described later in this paper, this accuracy is 
sufficient to observe the neutrino beam's bunch structure 
in the SK neutrino data. The system is described more 
completely in [33 |. 



IV. NEUTRINO INTERACTION SIMULATION 

The neutrino interaction simulation plays an impor- 
tant role both in estimating the expected number of neu- 
trino interactions and in deriving the energy spectrum of 
neutrinos from the data. The Monte Carlo program sim- 
ulates neutrino interactions with protons, oxygen, carbon 
and iron, which are the target materials of the neutrino 
detectors. 

In the simulation program, we include the follow- 
ing charged and neutral current neutrino interactions: 
quasi-elastic scattering {v N £ N'), single meson 
production {v N £ N' m), coherent tt production 
(f i6o(i2c,56Fe) £ TT i6o(i2c,56Fe)), and deep 

inelastic scattering {lyN ^ £ N' hadrons). In these 
reactions, N and N' are the nucleons (proton or neu- 
tron), £ is the lepton, and m is the meson. For the single 
meson production processes, the K and rj are simulated 
as well as the dominant tt production processes. If the 
neutrino interaction occurs in oxygen or other nuclei, the 
re-interactions of the resulting particles with the remain- 
ing nucleons in the nucleus are also simulated. 

A. Quasi-elastic scattering 

The formalism of quasi-elastic scattering off a free neu- 
tron used in the simulation programs is described by 
Llewellyn-Smith [35|. For scattering off nucleons in the 
nucleus, we use the relativistic Fermi gas model of Smith 
and Moniz [s^. The nucleons are treated as quasi- free 
particles and the Fermi motion of nucleons along with 
the Pauli exclusion principle is taken into account. The 
momentum distribution of the target nucleon is assumed 
to be flat up to a fixed Fermi surface momentum of 
225 MeV/c for carbon and oxygen and 250MeV/c for 
iron. The same Fermi momentum distribution is also 
used for all of the other nuclear interactions. The nu- 
clear potential is set to 27 MeV for carbon and oxygen 
and 32 MeV for iron. 



B. Single meson production 

Rein and Sehgal's model is used to simulate the reso- 
nance production of single tt, K and r] [s^, [H, [s^. This 
model divides the interaction into two parts. First there 
is the interaction 

I' + N £ + N*, 

which is then followed by 

N* -> 7r( or K or 77) + N', 

where N and N' are the nucleons, and N* is the baryon 
resonance like A(1232). The mass of the intermediate 
resonance is restricted to be less than 2 GeV/c^. To 
determine the direction of the pion in the final state. 
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we also use Rein and Sehgal's method for the dominant 
resonance P33(1232). For the other resonances, the di- 
rectional distribution of the generated pion is set to be 
isotropic in the resonance rest frame. The angular dis- 
tribution of has been measured for the vp — > ii^pn'^ 
mode [iOl and the results agree well with the Monte Carlo 
prediction. The Pauli blocking effect in the decay of the 
baryon resonance is taken into account by requiring that 
the momentum of the nucleon should be larger than the 
Fermi surface momentum. In addition, the delta may be 
absorbed by the nucleus. For these events there is no 
pion in the final state, and only a lepton and nucleon are 
emitted [4l|. We explicitly make this happen for 20% 
of the deltas produced. Single K and 77 productions are 
simulated using the same framework as for single tt pro- 
duction processes. 

Both the quasi-elastic and single-meson production 
models contain a phenomenological parameter (the axial 
vector mass. Ma), that must be determined by experi- 
ment. As the value of Ma increases, interactions with 
higher values (and therefore larger scattering angles) 
are enhanced. The Ma parameters in our Monte Carlo 
simulation program are set to be 1.1 GeV for both the 
quasi-elastic and single-meson production channels based 
on the analysis of the near detector data [3] . 

Coherent single tt production, the interaction between 
a neutrino and the entire nucleus, is simulated using the 
formalism developed by Rein and Sehgal [sst . Here, only 
the neutral current interactions are considered because 
the cross-section of the charged current coherent pion 
production was found to be very small at the K2K beam 
energy [i^ . 



Total o 



C. Deep inelastic scattering 



In order to calculate the cross-section for deep inelas- 
tic scattering, we use the GRV94 parton distribution 
functions [43! . Additionally, we have included the cor- 
rections in the small region developed by Bodek and 
Yang ■44]. In the calculation, the hadronic invariant 
mass, W, is required to be larger than 1.3 GeV/c^. Also, 
the multiplicity of pions is restricted to be larger than or 
equal to two for 1.3 < < 2.0 GeV/c^, because single 
pion production is already taken into account as previ- 
ously described. In order to generate events with multi- 
hadron final states, two models are used. For W be- 
tween 1.3 and 2.0 GeV/c^, a custom-made program [11] 
is employed while PYTHIA/JETSET ^ is used for the 
events whose W is larger than 2 GeV/c^. 

The total charged current cross sections including 
quasi-elastic scattering, single meson production and 
deep inelastic scattering are shown in Fig. [H] overlaid with 
data from several experiments. 
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FIG. 8: Charged current total cross section divided by E^, 
for neutrino nucleon charged current interactions. The solid 
line shows the calculated total cross section. The dashed, 
dot and dash-dotted hues show the calculated quasi-elastic, 
single-meson and deep-inelastic scatterings, respectively. 
The data points are taken from the following experiments: 
(A)ANL[4a|, (0)GGM77[4i], £.)GGM79(a) 49],(b) Mj, 
(*)Serpukhov51], (0)ANL82'52], (★)BNL86'4gj|, 

(■)CCFR90i53,], (T)CDHSW87^54^],^x)IHEP-JINR96t5|, 
(-|-)IHEP-ITEP79H, (□)CCFRR84[53|, and (a)BNL820. 



D. Nuclear effects 

The intra-nuclear interactions of the mesons and nu- 
cleons produced in neutrino interactions in the carbon, 
oxygen or iron nuclei are also important to consider for 
this analysis. Any absorption or change of kinematics 
of these particles will affect the event type classification. 
Therefore, the interactions of tt, K , rj and nucleons are 
also simulated in our program. These interactions are 
treated using a cascade model, and each of the particles 
is traced until it escapes from the nucleus. 

Among all the interactions of mesons and nucleons, 
the interactions of pions are most important, since both 
the cross sections for pion production for neutrino en- 
ergies above 1 GeV and also the interaction cross sec- 
tions of pions in the nucleus are large. In our simulation 
program, the following pion interactions in nucleus are 
considered: inelastic scattering, charge exchange and ab- 
sorption. The actual procedure to simulate these interac- 
tions is as follows: first the generated position of the pion 
in nucleus is set according to the Woods-Saxon nucleon 
density distribution 59). Then, the interaction mode is 
determined by using the calculated mean free path of 
each interaction. To calculate these mean free paths, we 
adopt the model described by Salcedo et al. [60|- The 
calculated mean free paths depend not only on the mo- 
mentum of the pion but also on the position of pion in 
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FIG. 9: The tt momentum distributions for neutral current 
single TT production processes oflt the water target and for 
the K2K neutrino beam at the near detector. The solid and 
dashed lines show the spectrum without and with the pion 
nuclear effects. 



V. THE FAR/NEAR FLUX RATIO 
A. Definition of the far-to-near ratio 

The effects of neutrino oscillation appear as a reduc- 
tion in the number of neutrino events and a distortion of 
the neutrino energy spectrum in SK. The observations for 
these quantities are compared to their expectations in SK 
to study neutrino oscillation. The ND measures the neu- 
trino flux and spectrum before neutrinos oscillate. Those 
measurements are then extrapolated by the expected ra- 
tio of muon neutrino fluxes at the far and near detector 
locations, the far-to-near (F/N) flux ratio, to predict the 
number of neutrino events and energy spectrum in SK. 

The neutrino flux at any distance from its source can 
be predicted when the geometry of the decay volume and 
the momenta and directions of the pion parents of neutri- 
nos are provided. Due to the finite size of the decay vol- 
ume and the detectors, the neutrino flux does not simply 
obey an rule (where L is distance from the neutrino 
source); rather the flux ratio between far and near detec- 
tors has some dependence on neutrino energy. Therefore, 
we define the F/N flux ratio, R^^^ , as 



the nucleus. 

If inelastic scattering or charge exchange occurs, the 
direction and momentum of pion are determined by us- 
ing the results of a phase shift analysis obtained from 
TT — N scattering experiments [6l| . When calculating the 
pion scattering amplitude, the Pauli blocking effect is also 
taken into account by requiring the nucleon momentum 
after the interaction to be larger than the Fermi surface 
momentum at the interaction point. 

This pion interaction simulation is tested by com- 
parison with data using the following three interac- 
tions: TT^^C scattering, tt^^O scattering and pion photo- 
production — > vr^ + X). The importance of in- 
cluding the proper treatment of nuclear effects is illus- 
trated in Fig. [HI which shows the momentum distribution 
for neutral current single vr" production in the water tar- 
get both with and without having them applied. 

The re-interactions of the recoil protons and neutrons 
produced in the neutrino interactions are also important, 
because the proton tracks are used to select quasi-elastic 
like events. This is done with the SciFi and SciBar 
near detectors, and allows us to estimate the neutrino 
energy. Nucleon-nucleon interactions modify the outgo- 
ing nucleon's momentum and direction, which also af- 
fects whether the nucleon will be above detection thresh- 
old [U . Both elastic scattering and pion production are 
considered. In order to simulate these interactions, a cas- 
cade model is again used and the generated particles in 
the nucleus are tracked using the same code as for the 



(6) 



where ^^^'-^^'> (E^) is the neutrino energy spectrum at 
SK (ND). 

The F/N flux ratio is estimated by our beam MC sim- 
ulation. In this simulation, while we use the Cho-CERN 
compilation as a reference model, we employ the HARP 
experiment [l6j result as an input for simulation of pion 
production. The pion production measurement done by 
HARP is of direct relevance for K2K, since it uses the 
same beam proton momentum and the same production 
target, and it covers a large fraction of the phase space 
contributing to the K2K neutrino flux. The details of 
the HARP measurements are described in Sec. IV Bl The 
pion monitor (PIMON) measurement is performed for 
a confirmation of the validity of the beam MC simula- 
tion. It gives us in-situ information on the momentum 
and the direction of pions entering the decay volume after 
they are focused by the horn magnetic fields although the 
PIMON is not sensitive to pions below 2 GeV/c (corre- 
sponding to neutrinos below 1 GeV) due to its threshold. 
A description of the PIMON measurement is given in 
Sec. Eel 



B. Prediction of far-to-near ratio from the HARP 
result 

The dominant uncertainty in neutrino flux predictions 
for conventional neutrino beams is due to the pion pro- 
duction uncertainty in the hadronic interactions of pri- 
mary beam protons with the nuclear target material. In 
this analysis, we use the results provided by the HARP 
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experiment at CERN as input to the pion production 
simulation. The HARP experiment precisely measured 
the positively-charged pion production in the interactions 
of 12.9 GeV/c protons in a thin aluminum target [iG^. 

The HARP experiment took data in 2001 and 2002 in 
the CERN PS T9 beamline, in order to study in a sys- 
tematic and accurate way hadron production for a vari- 
ety of produced hadrons (pions and kaons in particular) 
with large phase space coverage. Data were taken as a 
function of incident beam particle type (protons, pions), 
beam momentum (from 1.5 to 15 GeV/c), nuclear tar- 
get material (from hydrogen to lead), and nuclear target 
thickness (from 2% to more than 100% hadronic inter- 
action length fraction). Secondary tracks are efficiently 
reconstructed in the HARP forward spectrometer via a 
set of drift chambers located upstream and downstream 
with respect to a dipole magnet. Particle identification 
for forward tracks is obtained with a time-of-flight sys- 
tem, a Cherenkov threshold detector, and an electromag- 
netic calorimeter. 

In particular, the recent HARP pion production mea- 
surement [1(3] is directly relevant for the K2K F/N flux 
ratio because it is obtained for the same proton beam 
momentum (12.9 GeV/c) and nuclear target material 
(aluminum) as those used to produce the K2K neu- 
trino beam. Moreover, beam MC simulations show that 
the forward pion production region measured in HARP, 
30 < 6*^ < 210 mrad, 0.75 < < 6.5 GeV/c, matches 
well the pion production phase space responsible for the 
dominant fraction of the K2K muon neutrino fluxes at 
both the near and far detector locations. 

The result of the pion production measurements de- 
scribed in [l^ is incorporated into our beam MC simula- 
tion to estimate the neutrino spectra at ND and SK and 
the energy dependence of the F/N flux ratio in the ab- 
sence of neutrino oscillations. The relatively- normalized 
fluxes at ND and SK, <i)^^ and <i>^^, respectively, pre- 
dicted by HARP measurement, are shown in Fig. [TOl to- 
gether with the associated total systematic uncertainties, 
by the empty circles with error bars. Uncertainties in the 
primary and secondary hadronic interactions, in the pion 
focusing performance in the horn magnetic fields, and in 
the primary beam optics, are considered. Here, primary 
hadronic interactions are defined as hadronic interactions 
of protons with more than 10 GeV total energy in alu- 
minum, while secondary hadronic interactions are defined 
to be hadronic interactions that are not primary ones. In 
the following, the assumptions on systematic uncertain- 
ties affecting neutrino flux predictions are summarized. 

The uncertainty in the multiplicity and kinematics of 
TT"*" production in primary hadronic interactions is esti- 
mated based on the accurate HARP results. In this case, 
the HARP 7r+ Sanford-Wang parameters' uncertainties 
and correlations given in [1^] are propagated into flux 
uncertainties using standard error matrix propagation 
methods: the flux variation in each energy bin is esti- 
mated by varying a given Sanford-Wang parameter by a 
unit standard deviation in the beam MC simulation. An 
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FIG. 10: Relatively-normalized muon neutrino flux predic- 
tions at the near (top) and far (bottom) detectors. The empty 
circles with error bars show the central values and shape- 
only errors based on the HARP 7r+ production measurement, 
the empty squares with shaded error boxes show the central 
values and errors from the pion monitor (PIMON) measure- 
ment, and the dotted histograms show the central values from 
the Cho-CERN compilation of older (non-HARP) -k^ produc- 
tion data. The PIMON predictions are normalized such that 
the integrated fluxes above 1 GeV neutrino energy match the 
HARP ones, at both the near and far detectors. 



uncertainty of about 30% is assumed for the uncertainty 
in the proton-aluminum hadronic interaction length. The 
uncertainty in the overall charged and neutral kaon pro- 
duction normalization is assumed to be 50%. 

The systematic uncertainty due to our imperfect 
knowledge of secondary hadronic interactions, such as 
TT"*" absorption in the target and horns, is also consid- 
ered. We take the relatively large differences between the 
GCALOR/GFLUKA [UHl^ and GHEISHA [63^ de- 
scriptions of secondary interactions, also in comparison 
to available experimental data, to estimate this uncer- 
tainty. 
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FIG. 11: Prediction for the K2K muon neutrino F/N flux 
ratio in absence of osciUations. The empty circles with er- 
ror bars show the central values and systematic errors on the 
muon neutrino flux predictions from the HARP tv^ produc- 
tion measurement discussed in the text, the empty squares 
with shaded error boxes show the central values and errors 
from the pion monitor measurement, and the dotted his- 
tograms show the central values from the Cho-CERN com- 
pilation of older (non-HARP) tt"*" production data. 



interactions with the systematic error evaluation dis- 
cussed above, in the absence of neutrino oscillations, is 
shown in Fig. [11] as a function of neutrino energy. We 
estimate that the flux ratio uncertainty as a function of 
the neutrino energy binning used in this analysis is at 
the 2-3% level below 1 GeV neutrino energy, while it 
is of the order of 4-9% above 1 GeV. We find that the 
dominant contribution to the uncertainty in F/N comes 
from the HARP 7r+ measurement itself. In particular, 
the uncertainty in the flux ratio prediction integrated 
over all neutrino energies is 2.0%, where the contribu- 
tion of the HARP 7r+ production uncertainty is 1.4%. 
Table IIIII shows the contributions of all systematic un- 
certainty sources discussed above on the far-to-near flux 
ratio prediction for each neutrino energy bin. 

The dotted histograms in Figures [10] and [TT] show the 
central value predicted by using the "Cho-CERN" compi- 
lation for primary hadronic interactions, which was used 
in K2K prior to the availability of HARP data. In this 
case, the same Sanford-Wang functional form of pro- 
duction is employed to describe a CERN compilation of 
7r+ production measurements in proton-beryllium inter- 
actions, which is mostly based on Cho et al. data [15|]. A 
nuclear correction to account for the different pion pro- 
duction kinematics in different nuclear target materials is 
applied. The details of the Cho-CERN compilation are 
described in Sec. Ill CI We find that the predictions of 
F/N flux ratio by HARP and Cho-CERN are consistent 
with each other for all neutrino energies. Note that the 
difference between Cho-CERN and HARP central values 
represents a difference in hadron production treatment 
only. 



We account for the uncertainties in our knowledge of 
the magnetic held in the horn system. We assume a 
10% uncertainty in the absolute field strength, which is 
within the experimental uncertainty on the magnetic field 
strength and the horn current measured using inductive 
coils during horn testing phase Furthermore, a peri- 
odic perturbation in azimuth of up to ±15% amplitude 
with respect to the nominal field strength is assumed as 
the uncertainty in the field homogeneity, which is also 
based on the experimental accuracy achieved in the mea- 
surement of the mag netic field mapping in azimuth dur- 
ing horn testing |10 |. 

Finally, beam optics uncertainties are estimated based 
on measurements taken with two segmented plate ion- 
ization chambers (SPICs) located upstream of the tar- 
get. An uncertainty of 1.2 mm and 2.0 mrad in the mean 
transverse impact point on target and in the mean injec- 
tion angle, respectively, are assumed based on long-term 
beam stability studies [gj] . The uncertainty on the beam 
profile width at the target and angular divergence is also 
estimated, based on the ^--^20% accuracy with which the 
beam profile widths are measured at the SPIC detector 
locations [6^ . 

The F/N flux ratio, $SK/^ND^ predicted by the 
HARP 7r+ production measurement for primary hadronic 



C. Confirmation of far-to-near ratio by pion 
monitor measurement 

A conflrmation for the validity of the F/N ratio has 
been performed by in-situ pion monitor (PIMON) mea- 
surements. The PIMON was installed on two occasions 
just downstream the horn magnets to measure the mo- 
mentum (p^) versus angle (^^) 2-dimensional distribu- 
tion of pions entering the decay volume. The PIMON 
measurements were done twice: once measurement was 
done in June 1999 for the configuration of la period 
(200 kA horn current with 2 cm target diameter) and 
the other was done in November 1999 for the configura- 
tion of the other periods (250 kA horn current with 3 cm 
target diameter). 

A schematic view of PIMON is shown in Fig. (T^l PI- 
MON is a gas Cherenkov imaging detector which con- 
sists of a gas vessel, a spherical mirror, and an array of 
20 photomultiplier tubes. The Cherenkov photons emit- 
ted by pions passing through the gas vessel are reflected 
toward and focused onto the PMT array by the spher- 
ical mirror. Then, the PMT array on the focal plane 
detects the Cherenkov image. Due to the characteristics 
of the spherical mirror, photons propagating in the same 
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TABLE III: Contributions to the uncertainty in the far-to-near flux ratio prediction. The uncertainties are quoted in %. The 
six columns refer to different bins in neutrino energy, as shown in the table in units of GoV. 



Source 


0.0-0.5 


0.5-1.0 


1.0-1.5 


1.5-2.0 2.0- 


-2.5 


2.5- 


Hadron interactions 














Primary interaction rate 


0.3 


0.9 


0.9 


2.1 


0.2 


0.3 


tt'^ mult, and kinematics 


0.7 


2.0 


1.8 


2.1 


2.9 


4.7 


Kaon multiplicity 


0.1 


<0.1 


0.1 


<0.1 


0.1 


4.9 


Secondary interactions 


0.3 


1.2 


2.0 


2.1 


0.4 


0.7 


Horn magnetic field 














Field strength 


1.1 


0.8 


1.4 


4.2 


2.8 


3.9 


Field homogeneity 


0.3 


0.2 


0.5 


0.3 


0.6 


0.3 


Primary beam optics 














Beam centering 


0.1 


<0.1 


<0.1 


<0.1 


0.1 


0.1 


Beam aiming 


0.1 


<0.1 


<0.1 


0.1 


0.4 


0.2 


Beam spread 


0.1 


0.7 


1.7 


3.4 


1.0 


3.2 


Total 


1.4 


2.7 


3.6 


6.5 


4.2 


8.5 



To decay 
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FIG. 12: A schematic view of the pion monitor (PIMON). 
The PIMON consists of a gas vessel, a spherical mirror, and an 
array of 20 photomultiplier tubes. The gas vessel is filled with 
freon gas R-318 (G4F8). A wedge-shaped spherical mirror is 
set inside the gas vessel and Cherenkov light produced by the 
pions in the beam, represented by the concentric circles in the 
figure, is reflected by the mirror and directed to the array of 
photo-multiplier tubes which is set in the focal plane. 



direction are focused to the same position on the focal 
plane, giving us information on the direction of the pi- 
ons. The pion momentum is also obtained from the size 
of the Cherenkov ring. Furthermore, a momentum scan 
can be done by varying the refractive index of the inner 
gas. Therefore, the momentum and direction of pions 
can be measured separately by looking at the Cherenkov 
light distribution on the focal plane. 

As shown in Fig. I12i a wedge-shaped mirror is used as 
the spherical mirror to measure only 1 /30 of the beam 
assuming azimuthal symmetry of the distribution. Its 
top is aligned to be on the beam center. The reflection 



angle with respect to beam direction is 30°. 

An array of 20 PMTs (modifled R5600-01Q made by 
Hamamatsu Corporation) is set 3 m away from the beam 
center to avoid excess exposure to radiation. The size of 
the PMT outer socket is 15.5 mm in diameter and the 
sensitive area of the photocathode is 8 mm in diame- 
ter. They are arranged vertically at 35 mm intervals. 
The array can be moved by a half pitch of the inter- 
val along the array, and hence 40 data points (one point 
for every 1.75 cm) are taken for a Cherenkov light dis- 
tribution. The relative gain among 20 PMTs was cal- 
ibrated using Xe lamp before the measurements. The 
gain ratio between neighboring PMTs was also checked 
using Cherenkov photons during the run. The error on 
the relative gain calibration is estimated to be 10% for 
the June 1999 run and 5% for the November 1999 run. 
Saturation of the PMTs was observed in the June 1999 
run, which was corrected by a second order polynomial 
function. The uncertainty due to this correction was es- 
timated to be 4% [l3l- 

The gas vessel is filled with freon gas R-318 (C4F8). 
Its refractive index n is varied by changing the gas 
pressure using the external gas system. The data are 
taken at several refractive indices ranging between n = 
1.00024-1.00242 to make PIMON sensitive to different 
pion momenta. The refractive index was not adjusted 
beyond n = 1.00242 since the primary protons also emit 
Cherenkov photons when n exceeds this value, and be- 
come a severe background to the pion measurement. This 
corresponds to setting a momentum threshold of 2 GeV/c 
for pions, which corresponds to an energy threshold of 
1 GeV for neutrinos. The absolute refractive index is 
calibrated by the Cherenkov photon distribution from 
12 GeV primary protons with the refractive index set at 
71 = 1.00294. 

The Cherenkov light distribution for each refractive in- 
dex is taken by the PMT array. For the background sub- 
traction, a measurement with the mirror directed away 
from the direction of PMT array was performed. There 
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FIG. 13: Cherenkov light distributions for various refractive 
indices measured in November 1999. Dots show data and the 
histograms show the MC simulation. The refractive indices 
for each plot are as written in the figure. 



is still non-negligible background from electromagnetic 
showers which mainly come from the decay of neutral pi- 
ous, TT*' 27. The light distribution for this background 
is estimated using a MC simulation. The normalization 
in the subtraction is done by using the distribution mea- 
sured at the lowest refractive index, where the contribu- 
tion from the electromagnetic components is dominant. 
After all backgrounds are subtracted, the distribution of 
the Cherenkov light emitted from pions is obtained as 
shown in Fig. [131 The prediction of the MC simulation 
is superimposed as well. 

A x^-fitting is employed to extract the {p^^Ot^) 2- 
dimensional distribution from the Cherenkov light dis- 
tributions with various refractive indices. The (p^,07r)- 
plane is binned into 5 x 10 bins; 5 bins in above 
2 GeV/c with 1 GeV/c slice (the last bin is integrated 
over > 6 GeV/c) and 10 bins in from —50 mrad to 
50 mrad with 10 mrad slices. Templates of the Cherenkov 
light distributions emitted by pions in these bins are pro- 
duced for each refractive index using a MC simulation. 
Then, the weight of the contribution from each (pjr,^7r) 
bin being the fitting parameter, the MC templates are 
fit to observed Cherenkov light distributions. The fit- 
ting is done for the data in June 1999 and in November 
1999, separately. The resulting values of fitting parame- 
ters and errors on them in November 1999 run are shown 
in Fig. [H 

The neutrino energy spectra at ND and SK are de- 
rived by using the weighting factors obtained above and 
a MC simulation. The neutrino energy is binned into 6 
bins: 0.5 GeV bins up to 2.5 GeV, and integrated above 
2.5 GeV. The contribution of pions in each [p^^O^) bin 
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FIG. 14: The fitting results of pion (p^,^^) distribution in 
November 1999 run. The left figure shows the resulting cen- 
tral value of the weighting parameters, and the right figure 
shows the estimated fitting errors on them (no box means 
fitting errors are negligibly small.) 



to neutrino energy bins is estimated by a MC simulation, 
where to a good approximation it depends only on the 
pion kinematics and the geometry of the decay volume. 
Then, the neutrino spectrum is obtained by summing up 
these contributions weighted by fitted factors. Finally, 
the ratio of the neutrino spectra at SK to that at ND 
yields the F/N ratio. 

The extracted neutrino spectra and the F/N ratio from 
the PIMON data taken in November 1999 are shown 
in Fig. [TO] and [11] with empty squares and shaded er- 
ror boxes. All the systematic uncertainties in deriving 
them from the PIMON measurement are included in the 
errors, where the most dominant contributions to the er- 
ror on the F/N flux ratio come from the fitting error, the 
uncertainty in the analysis methodology, and the uncer- 
tainty in the azimuthal symmetry of the horn magnetic 
field. Further details on the systematic uncertainties in 
the PIMON measurement are described in h&i . 



D. The far-to-near ratio in K2K 

The F/N flux ratio used to extrapolate the measure- 
ments in ND to the expectation in SK is obtained in 
three independent ways: using the HARP measurement, 
the Cho-CERN model, and the PIMON measurement, as 
described in the previous sections. We find that all three 
predictions of the F/N ratio are consistent with each 
other within their measurement uncertainties. Among 
these measurements, we use the one predicted by the 
HARP measurement in our neutrino oscillation analysis 
described in this paper, since the HARP pion production 
measurement was done for the same conditions as K2K 
experiment: the proton beam momentum and the rele- 
vant phase space of pions responsible for the neutrinos 
in K2K are the same. In particular, the measured mo- 
mentum region by the HARP experiment reaches below 
2 GeV/c down to 0.75 GeV/c where the PIMON is insen- 
sitive. The HARP measurement also gives us the most 
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TABLE IV: Predictions for the F/N muon neutrino flux ratio 
as a function of neutrino energy, for the HARP model for tt^ 
production in primary hadronic interactions. The neutrino 
energy binning is also indicated. 



Energy Bin Number i 


[GeV] R^ 


(xlO-*^) 


1 


0.0-0.5 


1.204 


2 


0.5-1.0 


0.713 


3 


1.0-1.5 


0.665 


4 


1.5-2.0 


0.988 


5 


2.0-2.5 


1.515 


6 


2.5- 


1.720 



VI. MEASUREMENT OF NEUTRINO EVENT 
RATE AT THE NEAR DETECTOR 

The integrated flux of the neutrino beam folded with 
the neutrino interaction cross-section is determined by 
measuring the neutrino event rate at the near site. The 
event rate at the IKT detector is used as an input to the 
neutrino oscillation study. The stability of the neutrino 
beam is guaranteed by measuring the beam properties 
by the MRD detector. In addition, the LG and SciBar 
detectors measure the electron neutrino contamination in 
the beam to compare to our beam MC simulation. 



TABLE V; Fractional error matrix (SRiSRj) / (RiRj) obtained 
from the systematic uncertainties on the F/N flux predic- 
tions. The neutrino energy binning is the same as in Tab. II VI 
The values are given in units of 10"'^. 



Energy Bin 


1 


2 


3 


4 


5 


6 


1 


0.187 


0.002 


-0.036 


-0.372 


-0.281 


0.240 


2 


0.002 


0.728 


0.868 


1.329 


0.698 


-1.398 


3 


-0.036 


0.868 


1.304 


2.122 


1.041 


-2.040 


4 


-0.372 


1.329 


2.122 


4.256 


2.165 


-3.799 


5 


-0.281 


0.698 


1.041 


2.165 


1.779 


-2.678 


6 


0.240 


-1.398 


-2.040 


-3.799 


-2.678 


7.145 



accurate measurements on hadron production. 

The central values for the F/N flux ratio as a function 
of neutrino energy obtained from the HARP 7r+ produc- 
tion results, Ri, are given in Tab. IIVI where the index 
i denotes an energy bin number. The total systematic 
uncertainties on the F/N flux ratio as a function of neu- 
trino energy are given in Tab.fVl together with the uncer- 
tainty correlations among different energy bins, expressed 
in terms of the fractional error matrix (SRidRj) /{RiRj), 
where i,j label neutrino energy bins. The F/N central 
values and its error matrix are used in the analysis for 
neutrino oscillation described later. 

While the neutrino flux predictions given in this sec- 
tion are appropriate for most of the protons on target 
used in this analysis, a small fraction of the data was 
taken with a different beam configuration. The K2K-Ia 
period differed from the later configuration, as described 
in Sec. IIIBI As a result, the far/near flux ratio for June 
1999 is separately estimated, in the same manner as de- 
scribed above for later run periods. We find that the flux 
ratio predictions for the two beam configurations, inte- 
grated over all neutrino energies, differ by about 0.4%. 
The flux ratio prediction for the June 1999 beam con- 
figuration and the ND spectrum shape uncertainties are 
used to estimate the expected number of neutrino events 
in SK and its error for the June 1999 period. 



A. Neutrino event rate 

As described in Section III. A, the IKT water cherenkov 
detector serves to measure the absolute number of neu- 
trino interactions in the near site and to predict the num- 
ber of neutrino interactions in the far site. Since the 
IKT uses a water target and almost the same hardware 
and software as SK, the systematic error in the predicted 
number of interactions at the far site can be reduced. The 
intensity of the neutrino beam is high enough that multi- 
ple neutrino interactions per spill may occur in the IKT. 
When this happens it is difficult to reconstruct events. 
We employ Flash Analog- To-Digital Converters (FADCs) 
to record the PMTSUM signal (see Section III. A) and 
we can get the number of neutrino interactions by count- 
ing the number of peaks above a threshold. We set this 
threshold at 1000 photoelectrons (p.e.), approximately 
equivalent to a 100 MeV electron signal, to reject low 
energy background such as decay electrons from stopped 
muons. In Fig. 1151 the upper and lower figures show the 
number of peaks in a spill and the timing information 
of the peaks, respectively. We can clearly see the 9 mi- 
cro bunch structure of the beam in the lower figure. The 
fraction of multi-peak interactions in a spill is about 10% 
of single-peak spills. 

Sometimes the FADCs cannot identify multiple inter- 
actions if these events happen in the same bunch and 
the time gap between the interactions is too small. To 
correct for this possibility, we employ a MC simulation 
with multiple interactions to estimate the misidentifica- 
tion probability. The PMTSUM signals recorded by the 
FADC are simulated, so the same method can be used for 
the MC simulation and data. We found that the number 
of interactions in the fiducial volume is underestimated 
by 2.3% for multiple interactions. The multiple inter- 
actions contribute 34% of the total number of interac- 
tions, and we have to correct the number of events by this 
multi-interaction misidentification probability, according 
to Cmuiti = 1 + 0.023 X 0.34 = 1.008. 

The fiducial volume in the IKT is defined as a hor- 
izontal cylinder with axis along the beam direction (z- 
axis). The radius is 200 cm and the z coordinate is lim- 
ited to —200 cm < z < cm, where the center of the 
IKT ID is defined as z = cm, and the total fiducial 
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FIG. 15: The upper figure shows the number of neutrino inter- 
actions in a spiff. The fower figure shows the time distribution 
of peaks of PMTSUM signaf which are recorded by FADC. 
The beam's 9 micro bunch structure can be seen cfearfy. 



mass is 25 tons. The fiducial volume cut results in an 
almost pure neutrino sample, rejecting cosmic rays or 
muons generated by the beam in the materials surround- 
ing the IKT (beam-induced muons) . Figure [16] shows 
the vertex distributions of data and the MC simulation. 
Because we simulate only neutrino interactions without 
beam- induced muons and without the cosmic-ray muons, 
we can see excess events upstream of the z-distribution 
and the top part of the detector {y > 400 cm) in data. 
The data and the MC simulation are in good agreement 
in the fiducial volume. 

Two major background sources are considered. Cos- 
mic ray events usually have a vertex near the upper wall 
of the inner tank, but some events contaminate the fidu- 
cial volume due to failure of the vertex reconstruction. 
To estimate the background rate, we run the detector 
without the beam, replacing the spill trigger by a peri- 
odical clock signal. The beam-off data are analyzed in 
the same way as the neutrino data; it is found that cos- 
mic rays in the fiducial volume are 1.0% of the neutrino 
data. The other important background source is beam- 
induced muons which can be tagged by PMTs located in 
the outer detector. After the vertex cut, the remaining 
events are scanned with a visual event display and the 
fraction of beam induced muons found is 0.5%. In addi- 
tion, we had fake events which were produced by signal 
reflection due to an impedance mismatch of the cables 
in the 1999 runs only. The total background fraction is 
estimated to be 1.5% for runs starting in 2000, and 3.1% 
before 2000. 

The neutrino event selection efficiency is calculated 



FIG. 16: The upper figure shows the vertex distribution in z 
with an Rxy < 200 cm cut. The beam direction is defined as 
the z-axis. Crosses show data, boxes show the MC simulation 
normafized by area in the fiduciaf vofume, and verticaf bfack 
fines show the fiduciaf vofume. The middfe and fower figures 
show the X and y distributions, respectivefy. 



based on the MC simulation. The efficiency is defined 
as: 

iKT of reconstructed events in 25t) 

(# of generated events in 25t) 

Figure [T7] shows the selection efficiency as a function of 
neutrino energy. The overall efficiency including all ener- 
gies and all interaction types is 75% for the configuration 
with a 250kA horn current, and 71% for the 200kA config- 
uration. We had a problem with the FADC in November 
1999 which corresponds to 3% of all data, and the effi- 
ciency in this period was 5% lower than the other 250kA 
configuration periods. The dominant inefficiency comes 
from the single peak selection with a 1000 p.e. threshold 
by FADC. Figure [T5] shows the peak-finding efficiency of 
the FADC as a function of total charge. The 2000-2004 
data plotted in the figure shows that the efficiency curves 
are stable. The MC event selection efficiency is obtained 
by smearing the threshold assuming a gaussian distribu- 
tion, in which the mean and width are obtained by fitting 
the data. 

While the FADC can count the event multiplicity, they 
do not record information about each PMT channel. The 
ATM gives the timing and charge information of each in- 
dividual PMT channel which allows event reconstruction 
if the spill has a single interaction. Both the FADC and 
the ATM are required to derive the total number of neu- 
trino interactions: 
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FIG. 17: Neutrino selection efficiency as a lunction of neutrino 
energy (GeV). 
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FIG. 18: FADC peak finding efficiency is shown as a function 
of total charge (p.e.). Data taken in 2000, 2001, 2003, and 
2004 are shown with the MC simulation. 



^peak is the total number of PMTSUM signal peaks, 
is the number of single peak events, e^^""" is the 



single 



N' I, 

peak 

selection efheiency of neutrino events in the 25t fiducial 
volume, Rbkg is the background fraction, and Cmuiti is 
the multiple interaction correction. 

Table IVII shows the total number of neutrino interac- 
tions and the number of protons on target for each period. 



Table IVIII shows the systematic uncertainty on the 
number of neutrino interactions in the IKT. The domi- 
nant error is the uncertainty of the fiducial volume. From 
the comparison of neutrino interactions in data and the 
MC simulation, we quantitatively estimate the fiducial 
volume systematics. Varying the definition of the fidu- 
cial volume by about five times the vertex resolution, 
we observe a difference in the calculated event rate of 
1.8%. Most of the difference is due to the z-dependence 
of partially contained events. To estimate this systematic 
effect, we used deposited energy from neutrino interac- 
tions themselves for partially-contained events since the 
deposited energy is roughly linear in the distance from 
the vertex to the downstream wall of the ID. We use 
the "cosmic-ray pipe" muons, described in Section [ill A| 
to define the energy scale for partially-contained events 
within 2.3%. We do not see evidence for such a bias 
within the uncertainty of the energy scale. The fiducial 
uncertainty arising from a vertex bias is therefore 2.3%. 
We conservatively add those two numbers in quadrature 
to obtain 3.0% for the uncertainty of the fiducial volume. 

The energy scale uncertainty of the IKT is estimated 
by using cosmic ray muons which stopped inside of the 
detector and the reconstructed tt*^ mass which mostly 
comes from neutral current interactions. The absolute 
energy uncertainty of the IKT is estimated to be 
This energy uncertainty affects Nikt because of the 
FADC cut. We changed the threshold of the FADC and 
the effect of energy scale uncertainty is estimated to be 
0.3%. The FADC charge is calibrated by the total charge 
recorded by ATM using single interaction events with a 
lower threshold (200 p.e.). The stability of the charge 
scale of FADC is ±5% and its effect on Nikt is ±0.8%. 
Keak : ^Tat ^ud N^.Y depend on the FADC cut 
position, but N^^'^ should be independent of the cut if 
the efficiency correction is perfect. We calculate NjJ^'^ 
changing FADC cut from 200 to 2000 p.e. and confirm 
the total number of neutrino interaction is stable within 



TABLE VI: Number of neutrino interactions in the IKT. 



N, 



IKT 



= V, 



IKT 
obs 



ivrtotal 
peak 



single 



peak 



1 



clKT 



1 + Rbkg 



■ Cmuiti (8) 



where N}^"^ is the total number of neutrino interactions 
in the 25 t fiducial volume, N^^'^ is the number of events 
observed in the fiducial volume among single peak events. 



Period 




* obs 


7V'°*, 

'peak 


■\r single 
peak 


-'- ^ int 


la 


2.6 


4282 


109119 


89782 


7206 


lb 


39.8 


75973 


1854781 


1475799 


130856 


Ila 


21.6 


43538 


1061314 


832112 


73614 


lib 


17.1 


34258 


813599 


644723 


57308 


lie 


2.9 
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137533 


111834 


9346 
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B. Neutrino beam stability 
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FIG. 19: The upper figure shows the event rate in the IKT 
normahzed by the muon yield of the MUMON. One point 
corresponds to one month and the 2000 - 2004 data are shown 
in the figure. The middle line shows the average and the 
upper and lower lines show the RMS, 1.6%. The lower figure 
shows the IKT event rate divided by the MRD event rate, 
also showing the average and RMS (2.1%). 



TABLE VII: Systematic errors on iV, 
Source Error (%) 



IKT 

int 



Fiducial volume 


± 


3.0 


Energy scale 


± 


0.3 


FADC stability 


± 


0.8 


FADC cut position 


± 


1.5 


Event rate 


± 


2.0 


Background 


± 


0.5 


Multi- interaction 


± 


0.7 


Total 


± 


4.1 



±1.5%. In Fig. [191 the upper figure shows the IKT event 
rate normalized by muon yields at the SPD in the MU- 
MON (see Section II.A.2). We take 2.0% for the un- 
certainty of the event rate stability from the root-mean- 
square (RMS) of the distribution. A comparison between 
the event rates of the IKT and the MRD is shown in the 
lower figure of Fig. [19] as a consistency check. We assign 
statistical errors as systematic errors for the background 
and multiple interaction corrections due to the limited 
numbers of the sample. In total, we quote a ±4.1% error 
on the number of IKT neutrino events over the entire 
K2K run. 



The neutrino beam properties- the profile, beam di- 
rection, energy spectrum and event rate- are monitored 
by the MRD using interactions with iron in order to 
guarantee stability of the beam for the entire run period. 

A neutrino event is identified by a muon track in the 
MRD with the following selection criteria: only tracks 
within the time of the beam spill are accepted; tracks 
with a common vertex and a common timing are taken 
as a single event with the longest track assumed to be 
a muon; muons entering or exiting the detector are re- 
moved; and muons which traverse one or more iron plates 
> 1) are selected. 

The neutrino beam profile is obtained by measuring 
the vertex distribution in the MRD. The location of 
the center of the profile gives the beam direction. For 
this purpose, any muons with reconstructed energy lower 
than 0.5 GeV or higher than 2.5 GeV are rejected and 
a cubic fiducial volume of 6m x 6m in the upstream 9 
iron plates is used, which has a total mass of 419 tons. 
As shown in Fig. [201 the profile is well reproduced by the 
MC simulation. The profile center is plotted as a function 
of time in Fig. [^ The beam direction has been stable 
throughout the experiment within ±1 mrad from the SK 
direction. This satisfies the requirement for the direction 
to be within ±3 mrad, which is described in Section II. 
The stability of the profile width is also confirmed by the 
measurement. 

The muon energy and angular distributions are also 
continuously monitored. A cylindrical fiducial volume of 
radius 3m in the upstream 9 plates is used, where the 
mass is 329 tons. They show no change as a function 
of time beyond statistical uncertainty. The muon en- 
ergy distribution and the angular distribution are plot- 
ted monthly in Fig. [221 and Fig. [231 respectively. These 
results show that the t'^ energy spectrum is stable within 
2-4% depending on the energy bin. This is well within 
the spectrum error quoted in the spectrum analysis de- 
scribed in Section VII. 

In order to check the IKT event rate and to compare 
the neutrino cross sections in iron and water, the neu- 
trino event rate in the MRD iron is derived. The fiducial 
mass used is 72.8 tons, which is the 3m-radius cylinder 
of the upstream three iron plates only. For this purpose, 
Niayer > 2 is required in order to reduce the hadronic 
background. The stability of the IKT event rate di- 
vided by the MRD rate is plotted in the bottom figure of 
Fig.m 

Since the absolute normalization is less certain, 
a double ratio of the MRD and IKT event rates 
is calculated for the cross-section comparison be- 
tween iron and water. The result for this dou- 
ble ratio is (data/MC)MRD/(data/MC)iKT = 1-04 ± 
0.003(stat.) j:";"^(sys.). The average event rates for MRD 
and IKT data for the entire K2K run (except for K2K-Ia) 
are used. This double ratio should be unity if we correctly 
understand the iron/water neutrino cross-section ratio in 
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FIG. 20: Neutrino beam profiles measured by the MRD. 
Crosses show the measured profile and the boxes show the 
Monte Carlo prediction. Normalization is by number of en- 
tries. (Top left): horizontal profile for O.SGeV < < 
l.OGeV. (Top right): horizontal profile for l.OGeV < < 
2.5GeV. (Bottom left): vertical profile for O.SGeV < E^, < 
l.OGeV. (Bottom right): vertical profile for l.OGeV < E^, < 
2.5GeV. 
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FIG. 21: Stability of the neutrino beam direction measured 
by the MRD. The direction is plotted every five days for the 
entire experimental period. The solid line shows the SK di- 
rection and the dashed lines show ±1 mrad from the center. 
The direction is required to be stable within ±3mr. The top 
plot is for the horizontal direction; the bottom for the vertical 
direction. 



Muon energy spectrum stability 



TABLE VIII: Systematic errors of the event-rate double ratio 
(data/MC) MRD /(data/MC) IKT 



Source 


Error +(%) 


-(%) 


Fiducial volume 


+1.6 




-5.7 




Selection efficiency 


+1.2 




-5.7 




Tracking efficiency 


+1.0 




-1.0 




Beam direction 


+1.9 




-0.0 




MRD detector oriented total 


+2.9 


-8.1 


IKT detector oriented total 


+4.1 


-4.1 


Neutrino spectrum 


+0.9 




-0.9 




NC/CC ratio 


+4.0 




-3.7 




Non-QE/QE ratio 


+5.3 




-3.7 




Spectrum and neutrino int. 


+6.7 


-5.3 


Grand total 


+8.4 


-10.5 



the K2K energy range regardless of the absolute cross sec- 
tion. The sources of the systematic error are summarized 
in Table rvTITl Here, the uncertainties due to the neutrino 
spectrum, NC/CC ratio, and non-QE/QE cross-section 
ratio have the cancellation between IKT and MRD taken 
into account. 

The stability of the detector itself is confirmed as a 
whole by analyzing the off-spill data which is essentially 
cosmic-ray data taken between each beam spill. The 
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FIG. 22: Stability of the muon energy distribution measured 
by the MRD. Each bin is plotted every month for the entire 
experimental period except for K2K-Ia. The distributions are 
normalized by number of entries. 
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Muon angle distribution stability 
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FIG. 23: Stability of the muon angular distribution measured 
by MRD. Each bin is plotted every month for the entire ex- 
perimental period except for K2K-Ia. The distributions are 
normalized by their entries. 



event rate, angular distributions and Niayer 
are found to be stable. 



distribution 



C. Electron neutrino component 

The Vf, component in the beam is measured in the FGD 
system independently by the LG and by the SciBar de- 
tectors. In each detector we perform a search for Vg, in- 
teractions by looking for events with an electron in the 
final state. "Electron" events come essentially only from 
the Vf, component of the beam, since the Ve/v^i flux ratio 
is about 1.3 x 10~^, while the Ve/v^ flux ratio is about 
1.8 X 10~^, and the cross section for v^e scattering is 
about a factor 1.5x 10~^ smaller than that for CC scat- 
tering on a nucleon. The measurement of the v^. events 
validates the prediction of the Vf,jv^ flux ratio at the 
near location obtained from the beam MC simulation. 
The data-MC comparison also involves the cross-sections 
of and z/^ interactions, and the measurement is there- 
fore an important check of the MC simulation used to 
predict the number of interactions in SK. 

The LG measurement [65^ is performed by looking for 
i^e interactions taking place in the SciFi detector, with 
the electron detected in the downstream scintillator ho- 
doscope, and its energy measured in the LG calorime- 
ter. Electron events are selected by requiring the fol- 
lowing: an interaction vertex inside the SciFi fiducial 
volume, an energy deposit in the downstream scintilla- 



tor hodoscope system greater than 20 MeV (2.5 times 
larger than that expected from a muon), an energy de- 
posit in the LG greater than 1 GeV; and finally, no hit 
in the MRD matching the electron direction. For an ex- 
posure of 2.9 X 10^^ POT, 51 electron candidates are 
found with an estimated background of 24 induced 
events. The Ve/i'^i interaction ratio is estimated to be 
1.6 ± 0.4(stat.)lg'g(syst.)%, which is in agreement with 
the beam MC prediction of 1.3%. The dominant source 
of the t^e-induccd component of CC interactions is muon 
decay (87%) in the beam, and the remainder comes from 
kaon decay. 

The measurement of the Vf, contamination in the beam 
has also been performed using the SciBar plus EC detec- 
tor, with a statistics corresponding to 2.1 x 10^^ POT. 
The search for electrons is mainly based on the signals 
from the EC, the electromagnetic calorimeter which fol- 
lows the tracking section of the SciBar. The two planes of 
the EC correspond to approximately 5.5 radiation lengths 
(Xq) each, so that electrons in the 1 GeV energy range 
are almost fully contained. The average energy lost by 
a muon or pion traversing one plane is small, of the or- 
der of 50 MeV (with - 60 MeV fuU width). Therefore, 
to look for electrons we select events with a large signal 
in a restricted region of the first EC plane. We require 
El > 350 MeV, where Ei is the energy in a cluster of 
20 cm width centered in the module with maximum sig- 
nal. We also apply conditions on the energy release E2 
in the 20 cm wide cluster of highest energy of the second 
plane of the EC. We require E2/E1 to lie in the 0.2-1.1 
interval, which from MC simulation we know to contain 
95% of the electron events. In the tracking volume of the 
SciBar we search for a reconstructed track pointing to 
the selected high energy clusters in the two projections 
of the EC. The efficiency for reconstructing the electron 
as a track in the SciBar is high, given the low density of 
the detector (Xq = 40 cm). Finally, we impose a fidu- 
cial volume cut on the interaction vertex, defined as the 
starting point of the electron track, and we also require 
that outside the selected cluster the energy in each EC 
plane does not exceed 30 MeV. 

The selected sample consists of 42 electron candidate 
events. The visual examination of the display of these 
events allows us to discard 9 events, easily identified as 
background from neutrino interactions originating out- 
side the fiducial volume of the SciBar, or as interactions 
with 7r° production. From the number of ir^ events iden- 
tified at the scanning level we can estimate, with a cor- 
rection obtained from a MC simulation, the number of 
tt'^ events which cannot be distinguished from electron 
events, and constitute an irreducible background. Our fi- 
nal sample, with 90% electron purity, contains 33 events, 
with a background estimate of 3 ± 2 events. The charac- 
teristics of the selected events are compared to those of 
MC events resulting from a full simulation of Ve interac- 
tions in the SciBar and including also a 10% background 
from interactions. The data-MC comparison is fair, 
except for an excess of high energy electrons in the data. 
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FIG. 24: Electron energy spectrum for the candidates of in- 
teractions selected in SciBar. The MC spectrum (histogram) 
is normalized to the data separately for signal (30 events) and 
background (3 events). The contents of the highest energy bin 
is integrated between 3.6 and 8 GeV. 



distributions of versus 9^ are used to measure the neu- 
trino energy spectrum. The spectrum is fitted by using 
a method to compare observed (p^, distributions 
to MC expectations. 



A. The fitting method 

In order to obtain the neutrino energy spectrum, the 
(Pa" ^m) distribution is fit with the MC expectation as 
shown in Fig. 1251 The neutrino energies are divided into 
eight bins as shown in Table IIXI For the MC expec- 
tation, the (j)^,9^) distribution is prepared for each E,y 
bin and separately for QE and nonQE interactions; 8x2 
distributions are prepared in total for each event sample. 

The free parameters in the fit arc the neutrino energy 
spectrum parameters for eight energy bins {ff,...,f^) 
and a parameter, i?nQE, which represents the relative 
weighting of CC-nonQE events to CC-QE events. The 
systematic uncertainties, such as nuclear effects, the en- 
ergy scale, the track finding efficiency, and other detector 
related systematics, are also incorporated as the fitting 
parameters (/). The contents in (TO,n)-th bin of the 
(Pm' ^m) distribution, N^^, is expressed with the 16 tem- 
plates and the fitting parameters as, 



The energy spectrum of the electrons is shown in Fig. [24l 
The electron energy is obtained by correcting the energy 
measured in the EC for the energy lost in the tracking 
section of the SciBar and for the longitudinal leakage, the 
average correction being of the order of 20%. The excess 
may indicate an underestimate of kaon production in the 
beam simulation, but the statistics is too small to draw 
firm conclusions. Finally, we use the MC simulation to 
extrapolate our measurement, which is only sensitive to 
electron energies larger than 500 MeV, to the full energy 
range. Our result for the interaction ratio v^CC / u^^CC 
is 1.6 ± 0.3(stat.) ± 0.2(syst.)%, consistent with the MC 
prediction of 1.3%. 

The consistency of the measurements in the LG and 
SciBar between themselves and with the MC predictions 
confirms the quality of the measurements and of the MC 
simulation. However, since the measurements are lim- 
ited to a restricted energy region, for all the analyses 
described in the paper we use the Ve component as given 
by the MC simulation. 
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where P, N^'^^f^^ and N^^[^'^^^ are a normaliza- 
tion parameter, the number of expected contents in the 
(to, n)-th bin for QE interaction and that for non-QE in- 
teraction for the i-th neutrino energy bin. We take the 
between the observed distributions, N^^^, and iV^^. 
During the fit, the flux in each energy bin and -RnQE 
are re-weighted relative to the nominal values in the MC 
simulation. The flux for Ei, = 1.0 — 1.5 GeV bin is fixed 
to unity for the normalization, and another set of pa- 
rameters is prepared for relative normalization of each 
detector. 

The functions are separately defined for each detec- 
tor and then summed to build a combined function 
as 



VII. MEASUREMENT OF NEUTRINO 
SPECTRUM AT THE NEAR DETECTOR 



The neutrino energy spectrum before oscillation is 
measured with near detector IKT, SciFi and SciBar CC 
event samples. The neutrino energy is reconstructed from 
the muon kinematics parameters and 0^ assuming a 
QE interaction as given in Eq. [51 The two-dimensional 



2 _ 2 I 2 , 2 
XnD — XlKT + XSF + XsB 



(10) 



Finally, a set of the fitting parameters (/,f,i?nQE : /) is 
found by minimizing the function. The best fit values, 
their error sizes and the correlations between them are 
used as inputs to the oscillation analysis, as described 
in Sec IIX Al The following subsections will describe the 
definition of x^ for each subdetector and the results of 
fit. 



25 



TABLE IX: The E,, interval of each bin. 





r(p rip r(p rip rip rp rip rip 
Jl J2 JZ Ji J5 J(i J7 Ja 


[GeV] 


0.0-0.5 0.5-0.75 0.75-1.0 1.0-1.5 1.5-2.0 2.0-2.5 2.5-3.0 3.0- 



TABLE X: The summary table for the number of observed 
events in the FClR/i sample, the efficiency and the purity of 
the CCQE events estimated with the MC simulation. 

# of events CCQE efficiency (%) CCQE purity (%) 
FClR^i 52TTO 53?7 57^9 



B. Definition of x for IKT 

To measure the neutrino energy spectrum, we select a 
QE enriched data sample called the fully contained one- 
ring ^-like (FCIR^) sample. In addition to the IKT 
event rate selection (see Section VLA), we require four 
additional conditions. These are that all visible parti- 
cles are inside the IKT detector (fully contained), one 
Cherenkov ring is found (one-ring), the particle identi- 
fication is a muon (/x-like), and the reconstructed muon 
momentum is greater than 200 MeV/c. This last condi- 
tion is to ensure the quality of the event reconstruction. 
After these cuts, the fraction of CCQE events is about 
60%. See Table 1X1 for the data summary and CCQE frac- 
tion. The requirement of full containment in the IKT 
suppresses events at high momentum compared to the 
other detectors. Figs. [26l [27l and [28] show the ring num- 
ber likelihood, particle identification likelihood, and the 
fully contained versus partially contained event (FC/PC) 
separation, respectively, used in the FClR/x event selec- 
tion. Any discrepancies between data and the MC simu- 
lation observed are used to estimate some of the system- 
atic errors in the IKT mentioned below. 

The selected data and the MC simulation are then 
binned into 2-dimensional distributions of muon momen- 
tum versus the scattering angle (0^^). The momentum is 
divided into 16 100 MeV/c bins from - 1600 MeV/c. 
The scattering angle is divided into 10 bins where the 
first nine are in increments of 10 degrees from 0°-90° and 
the final bin contains all events with the angle greater 
than 90°. The bins for which the expected number of 
events is greater than 10 are used for the analysis (80 
bins in total). 

The neutrino spectrum is derived by comparing the 
observed data and the weighted sum of Monte Carlo ex- 
pectations using a test. The is defined as: 
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where N^'^ is the number of observed events for data 
for (m, n)-th bin, N^'^ is the number of MC events as 
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FIG. 25: Schematic view of the binning of the data and Monte 
Carlo events for the spectrum fit. The left plot shows a corre- 
lation between and for FClR/i events in the IKT data 
used for the spectrum fit. The right plots show those for the 
MC sample separately prepared for each neutrino energy bin 
and for QE and non-QE interactions. 



given by equation [SI e is the fitting parameter for en- 
ergy scale (e = 1 is a nominal value) which scales the 
muon momentum, am.n is the error including statistical 



and systematic errors, and a, 



energy 



is the estimated un- 



certainty of the energy scale, -1-3/— 4%. 

The parameters ff, i?nQE, P, and e are the fitting pa- 
rameters as discussed above. 

The systematic errors that are introduced in the fit 
come mainly from the event selection. Those errors are 
the ring counting likelihood, the particle identification 
likelihood, the FC/PC cut, the event vertex and direc- 
tion reconstruction, and the energy scale. The other sys- 
tematic errors are from the detector calibration and the 
axial vector mass (Ma), which is used for our neutrino 
interaction model. 



C. Definition of for SciFi 

The SciFi tracking detector can observe charged parti- 
cle tracks produced in neutrino- water interactions. Since 
SciFi uses the same target material (water) as IKT and 
SK, systematic uncertainty due to different target nuclei 
is reduced. The data taken during K2K-Ib and K2K-IIa 
periods have been analyzed. 

We choose charged current events in which a muon 
track starts from an interaction in the SciFi fiducial vol- 
ume and stops in the MRD detector. The fiducial region 
is defined as a rectangle 1.1 m to each side of the detec- 
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FIG. 26: The distribution of ring counting likelihood for the 
IKT. Those events that have a likehhood less than or equal 
to 0.0 are considered to have one ring; those above 0.0 are 
considered to be multi-ring. In this plot, data are the circles 
and the MC simulation is the histogram. The hatched his- 
togram shows the CCQE component. Only statistical errors 
are shown for data. 



10000 



5000 




PID iilcelihood 

FIG. 27: The distribution of particle identification likelihood 
for the IKT. The events with a likelihood greater than 0.0 are 
fi-like while those less than or equal to 0.0 are e-like. Data 
are the circles and the MC simulation is the histogram with 
the CCQE component shown as the hatched area. 



tor's center in both x and y, covering the first to 17*'' 
water containers. The fiducial mass is 5.59 ± 0.07 metric 
tons. The primary (muon) track should match a hit in 
the downstream scintillator hodoscope, with no matching 
hit in the upstream one. It should also match a recon- 
structed track which has at least penetrated one piece of 
steel and produced hits in two layers in the MRD. The 



FIG. 28: The largest charge in a PMT for a IKT event. 
Events that have a charge less than 200 p.e. are considered FC 
events. The rest are PC events. Data are the circles and the 
MC simulation is the histogram, with the CCQE component 
shown as the hatched area. 





lb 


Ila 


Detector Configuration 


SciFi-fLG-fMRD 


SciFi-fMRD 


POTsciPi (xlO^**) 


39.70 


22.04 


SciFi-MRD (track match) 


6,935 


5,188 


SciFi-MRD(cell match only) 


1,403 


1,743* 


SciFi-LG 


2,666 





TABLE XI: Number of events in each event category observed 
in SciFi. *The events are not used for the spectrum recon- 
struction. 



muon momentum threshold is 675 MeV/c for period-lb 
and 550 MeV/c for period Ila. 

For the K2K-Ib data, to enhance sensitivity to low 
energy neutrino interactions, we also include events in 
which the muon stops in the Lead Glass calorimeter (LG) 
and events which penetrate even single active layer in 
the MRD. Here, the response of the matched LG cluster 
should be greater than 100 MeV, which prevents proton 
tracks from being identified as muon tracks. By includ- 
ing these samples the threshold for period-lb data is 
reduced to 400 MeV/c. 

The number of events for each event category, together 
with the total POT of the beam spills used for the analy- 
ses, are summarized in Table IxH In total 17,935 charged 
current candidate events have been collected for the data 
corresponding to POTsciFi=6.17x 10^^ POT. After cor- 
recting for the changing detector configuration and effi- 
ciency, the event rate per POT was stable from month 
to month. 

In the spectrum reconstruction, there must be one or 
two tracks in the event, including the muon; events with 
three or more tracks are discarded, which amounts to 
about 3% of the total. The two-track sample is further 
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TABLE XII: QE efficiency and QE purity of all events for 
K2K-I and K2K-IIa for each SciFi subsample, estimated from 
the MC before fitting. 





QE efficiency (%) 


QE purity (%) 




K2K-I K2K-IIa 


K2K-I K2K-IIa 


1-track 


39 36 


50 57 


QE 


5 5 


53 58 


nonQE 


2 2 


11 12 


total 


46 42 





subdivided into a QE and a non quasi-elastic sample. For 
QE events, the muon momentum and muon angle are suf- 
ficient to predict the angle of the recoil proton because of 
the simple two-body kinematics. If the observed second 
track matches this prediction, within 25 degrees, then it 
is likely the event is QE. If the second track is more than 
25 degrees from the predicted angle, then it is very likely 
not a QE event, as shown in Fig. [29l 
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FIG. 29: Example of the distribution of /S.6p, the difference 
between the observed and predicted (assuming QE interac- 
tion) angle of the second track. This distribution is for SciFi 
2-track events for the K2K-I period. The data (circles) and 
MC (histogram) are compared, and the shaded region shows 
the QE fraction in the MC. 

The above data sample gives three subsamples each for 
K2K-I and K2K-IIa. Tab. Kn] shows the QE efficiency 
(fraction of QE MC events ending up in each subsample) 
and QE purity (fraction of events in the subsample which 
are QE) for K2K-I and K2K-IIa. The LG stopping events 
are separated from the rest of the K2K-I data, giving an- 
other three subsamples - a total of nine. These samples 
are then divided into seven angle bins from to 60 de- 
grees, in 10 degree increments, and also into eight muon 
momentum bins, which have the same intervals as the 
neutrino energy bins in Tab. lIXI The kinematics of these 



interactions are such that the bins at both high angle and 
high momentum contain neither data nor MC events and 
are not used in the analysis. The high muon thresholds 
also mean there are no data in the lowest energy region, 
and it is also not used in the fit. 

In addition to the energy spectrum and RnQB param- 
eters, the fit includes several systematic parameters spe- 
cific to the SciFi detector. Three refer to uncertainties in 
the reconstructed energy of each event. Since the muon 
momentum is reconstructed from the range of the muon, 
uncertainty in the material assay, and in the dE/dx used 
in the GEANT Monte Carlo simulation, as well as the un- 
certainty in the outgoing muon energy from a neutrino 
interaction, are important. A muon energy scale param- 
eter P|^aio '^ith an uncertainty of ±2.7% is applied to 
the measured muon momenta. The second parameter 
-^LG-donsity accouuts for an uncertainty of ± 5% in the 
muon energy loss specifically when it passes through the 
LG, and we consider this to be an uncertainty in the LG 
density. A final energy scale applies only to the energy 
reconstructed from the visible energy clusters in the LG 
for LG stopping events. This is expressed as the param- 
eter -pLG-ciustor ^'^d is given as an energy shift in GeV. 
The uncertainty used in the fit for the LG cluster energy 
calibration is ± 30 MeV. The muon momentum for LG 
stopping events is reconstructed only from the range in 
SciFi and the cluster energy observed in the LG. 

There are two other systematic parameters: a 
migration between one-track and two-track events, 
-f:fnd-track-off' which accounts for the tracking efficiency 
for short, second tracks; and a migration between the 
two-track QE and non-QE samples, -Prcscattoringi ac- 
count for nuclear effects such as the uncertainty in the 
re-scattering of the proton as it leaves the nucleus from 
a QE interaction. Both parameters are simple linear mi- 
grations, and we take a 5% uncertainty in the former and 
a 2.5% uncertainty in the latter. They move a fraction 
of events from a particular and 0^ bin from one of the 
three subsamples to the same bin in another of the three 
subsamples. 

There are 286 bins of SciFi data included in the fit, 
and six SciFi specific systematic parameters including 
normalization ^Noim- C)ur is defined as the negative 
of the logarithm of the Poisson likelihood for the binned 
data, plus terms arising from the pull of five system- 
atic errors (but not normalization): 

= -21nA(0) 

= 2 _ ^obs^ ^ ^obs^ ln«y<?,(0))] 

XEscalc ^ XlG density XlG cluster energy 
+ X2nd track eff + XRescattoring- (12) 

in which N^^{0) and N^^"^ are the predicted and ob- 
served values in the (77i,n)-th bin for some values of the 
parameters 9. This is the simplified version given in 
Ref. [2l| and its minimum follows a distribution. To 
this, we add terms arising from the systematic errors. 
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TABLE XIII; Number of observed events in each event cate- 
gory and the efficiency and purity of the QE events estimated 
with the MC simulation for the SciBar detector. 





# events 


QE efficiency (%) 


QE purity (%) 


1-track 


7256 


50.0 


57.8 


QE 


1760 


15.4 


71.3 


nQE 


2014 


3.7 


15.9 


Total 


11030 


69.1 





a Gaussian distribution. The Xdist ^^^^ defined as: 
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D. Definition of for SciBar 



The normalization parameter of each sample is given 



by 



The CC event selection in the SciBar detector is simi- 
lar to that of the SciFi detector. The fiducial volume of 
SciBar is defined as a rectangle that extends 2.6 x 2.6 yc? 
around the beam axis, from the second to 53'"'' layer 
of scintillator. The fiducial mass is 9.38 tons. Events 
with any track starting in the SciBar fiducial volume 
and matched with a track or hits in the MRD are se- 
lected as CC candidates (SciBar-MRD sample). This 
requirement imposes a threshold of 450 MeV/c on p^, 
reconstructed from its range through SciBar and MRD. 
According to the MC simulation, 98% of the events se- 
lected by this requirement are CC induced events, and 
the rest are neutral current (NC) interactions accompa- 
nied by a charged pion or proton which penetrates into 
the MRD. The resolution, scale uncertainty and 0^ 
resolution are 80 MeV/c, 2.7% and 1.6 degrees, respec- 
tively. 

Events with one or two reconstructed tracks are se- 
lected from the SciBar-MRD sample. The two-track 
events are subdivided into two categories — QE samples 
and nonQE samples — by using kinematic information, 
A0p, like SciFi. The number of observed events and the 
fraction and efficiency of QE interactions estimated with 
the MC simulation for each event category are summa- 
rized in Table Kml 

The SciBar term of the consists of two components, 
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rjNorm 
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pSB 

^2trk/ltrk 
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pNorm 
^nonQE 


_ pSB 
~ ^Norm 


pSB pSB 
^2trk/ltrk ' -^nonQE/QE 


(18) 



where ^Norm the overall normalization factor, and 
^'2^t?k/itrk and iff„QE/QE are the parameters to vary 
the ratio of the number of 2-track events to that of 1- 
track events and the ratio of the number of CC-nonQE 
events to that of CC-QE events, respectively, within their 
systematic uncertainties. We consider the uncertainties 
from nuclear effects and detector systematics. The possi- 
ble difference between carbon and oxygen is included in 
the nuclear effect uncertainty. Because the nuclear effects 
are common source of the uncertainties for f'2tri{/itrk and 
P^^„„,„^, their correlation is also estimated. In addi- 

tion, -Pp^scaio introduced to account for the uncertainty 
of the energy scale of the muon reconstruction. The mo- 
mentum scale uncertainty is 2.7% as described above. 

The Xsyst is calculated with constraint parameters, in- 
cluding their correlation: 



Xsyst 



(19) 



where -Pgyst represents a set of systematic parameters, 
Pq are their nominal values, and V is a covariance ma- 
trix. Three systematic parameters , Pp^^gcaie' -f2^t?k/itrk' 
and ^n(^iQE/QE are included in Pgyst; they are defined 
as relative weighting factors to the nominal MC expec- 
tation and all components of Pq are set to unity. The 
uncertainties and correlation among the parameters are 
evaluated to be 



The Xdist is calculated by the binned likelihood method 
using the (p^,0i^) distribution. The bin widths are 
0.1 GeV/c for and 10 degrees for B^. Bins with the 
expected number of events greater than five are used for 
the fit. In total, 239 bins are used for the analysis. The 
number of observed events in each bin, A^^*^^, is assumed 
to follow a Poisson distribution with the mean N^^. A 
bin-by-bin systematic uncertainty on hadron contamina- 
tion is implemented by the convolution of a Poisson with 
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+ (0.059)2 +(0.017)2 
+ (0.017)2 +(0.058)2 

The dominant error sources are track finding efficiency 
for -P|t^i^/itrk and the nuclear effects uncertainty for 

pSB 

-^nonQE/QE- 
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E. Fit results 

The minimum ^ point in the multi-parameter space 
is found by changing the spectrum shape, i?nQE and the 
systematic parameters, where the MINUIT program h- 
brary j66j is employed. The central values and the errors 
of the fitting parameters are summarized in Table IXIVI 
All the systematic parameters stay within their estimated 
errors. The result of the spectrum measurement is shown 
in Fig.[3niwith the prediction of the beam MC simulation. 

The results of the measurements with individual de- 
tector data are also shown in Table IXTVl In the fit with 
only IKT data, the energy spectrum parameters are fixed 
to their default values for the high energy region where 
there is little or no acceptance. For the same reason, 
the low energy region is fixed for SciFi and SciBar. All 
the fitting parameters are in good agreement within their 
errors each other except for i?nQE- 

The p^, Q^i and g^ec distributions for the IKT, SciFi 
and SciBar samples are shown in Fig. [^ [551 In these 
figures, the reconstructed Q"^ distributions((7j^(,j.) are con- 
structed by assuming that the interaction was CC-QE 
and using the reconstructed energy under this assump- 
tion. The expected distributions of the MC simulation 
with the best-fit parameters are also shown. The MC 
simulation reproduces all the distributions well. 

The discrepancy in i?nQE is treated as a systematic er- 
ror. However, the value of i?nQE is strongly correlated 
with the Ey spectrum as well as the other systematic pa- 
rameters such as ^'nc^QE/QE' Order to evaluate i?nQE 
with each detector data set under identical fitting con- 
dition, a second fit is performed. In the second fit, the 
Ey spectrum and the systematic parameters, except for 
the overall normalization, are fixed at the best fit values 
obtained with all the three detectors. The best fit value 
of -RnQE for each detector in the second fit is (IKT, SciFi, 
SciBar) = (0.76, 0.99, 1.06), respectively, while the fit re- 
sult with three detectors is 0.96. Therefore, an additional 
error of 0.20 is assigned to i?nQE in order to account for 
the discrepancy. 

The errors of the measurement are provided in the form 
of an error matrix. Correlations between the parameters 
are taken into account in the oscillation analysis with this 
matrix. The full elements in the error matrix are shown 
in Table KVl 



VIII. SK DATA 

A. SK data reduction 

Beam-induced neutrino events in SK are selected ac- 
cording to criteria described in this section. Following 
selection, events which are fully contained in the SK fidu- 
cial volume are reconstructed using similar methods to 
the SK atmospheric neutrino analysis and then used in 
the K2K oscillation analysis. 
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FIG. 30; The neutrino energy spectrum measured at the near 
site, assuming CC-QE. The expectation with the MC simula- 
tion without reweighting is also shown. 



In order to select those neutrino interactions which 
come from the accelerator at KEK, two Universal Time 
Coordinated time stamps from the worldwide GPS sys- 
tem are compared. Both Tkek for the KEK-PS beam 
spill start time, and Tqk for the SK trigger time are 
recorded. The time difference AT = Tsk~Tkek~TOF , 
where TOF is a time of fiight, is distributed from and 
l.I /isec matching the timing width of the beam spill of 
the KEK-PS. The maximum difference of the synchro- 
nization between two sites is measured to be less than 
200 ns by using an external atomic clock. For this reason 
we require the AT for selected events to be between -0.2 
to 1.3 ^ sec. 

In addition to the timing criteria, the following cuts 
are required: 

1. In order to remove decay-electrons from the sample, 
events must have no activity in the 30 /iS before the 
event. 

2. There must be a minimum number of photo- 
electrons seen within a 300 ns timing window. 
The required number of photo-electrons are 200 for 
K2K-I and 94 for K2K-II. 

3. Fully contained events are selected by requiring no 
activity in the outer detector. 

4. A selection is made to remove events with PMTs 
which sometimes begin to produce light because 
of a discharge around the dynode. These events 
have easily identified characteristics such as a tim- 
ing distribution which is much broader than neu- 
trino events, and a repeating pattern of light in the 
detector. 

5. At least 30 MeV energy must be deposited in the 
inner detector. 
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FIG. 31: The distributions for each event sample of all near detectors with the MC simulation after fitting, given by open 
histograms. The hatched areas are the CCQE components in the MC distributions. 
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FIG. 32: The 9^ distributions for each event sample of all near detectors with the MC simulation after fitting, given by open 
histograms. The hatched areas are the CCQE components in the MC distributions. 
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FIG. 33: The g^cc distributions for each event sample of all near detectors with the MC simulation after fitting, given by open 
histograms. The hatched areas are the CCQE components in the MC distributions. 
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TABLE XIV: Results of the spectrum measurement. The best fit value of each parameter is listed for the fits with all the 
detectors' data, with the IKT data, with the SciFi data and with the SciBar data, respectively. The reduced (Xtotai/DOF) 
and the averaged x'^ of each detector (x'^/Nbin) are also shown. 



parauictcr 


Combined 


iKT only SciFi (nily SciBar only 


/i (0.00-0.50 GeV) 
/2 (0.50-0.75 GeV) 
js [\j . i 0- 1 .1)1) Kjrew ) 
/4 (1.00-1.50 GcV) 
/s (1.50-2.00 GcV) 
/e (2.00-2.50 GeV) 
/t (2.50-3.00 GeV) 
fs (3.00- GeV) 

RnQE 


1.657 ±0.437 
1.107 ±0.075 

1.104I: It U.UDl 

= 1 

0.911 ±0.044 
1.069 ± 0.059 
1.152 ±0.142 
1.260 ±0.184 
0.964 ±0.035 


2.372 ± 0.383 =1 =1 
1.169 ±0.072 0.882 ±0.317 1.166 ±0.251 

i.UOi zt U.UDD l.±Dl ± U.ZUi l.±4D ziz U.io4 

=1 =1 =1 

0.709 ±0.151 0.980 ±0.107 0.963 ± 0.070 
= 1 1.188 ±0.096 0.985 ±0.086 
= 1 1.062 ±0.230 1.291 ±0.283 
= 1 1.323 ±0.203 1.606 ±0.749 

0.589 ±0.071 1.069 ±0.060 1.194 ±0.092 


plkt 

^ Norm 
plkt 
^ energy 


0.948 ± 0.024 

0.984 ±0.004 


1.172 ± 0.046 — — 
0.993 ± 0.007 — — 


pHF 

Norm 
pSF 

pSF 

LG — density 

PLG-cluster [GeV] 
pSF 

'^2nd-track-eff 
pSF 

rescatterinff 


1.009 ±0.029 
0.980 ± 0.006 
0.929 ±0.012 
-0.001 ± 0.002 
0.959 ±0.014 
1.048 ± 0.055 


— 0.925 ±0.058 — 

— 0.980 ±0.007 — 

— 0.928 ±0.012 — 

— -0.002 ± 0.003 — 

— 0.932 ±0.017 — 

— 0.993 ± 0.062 — 


Norm 
pSB 

p — scale 
pSB 

-'^2trk/ltrk 
pSB 

-f^noutJE/CJE 


0.998 ± 0.010 
0.976 ± 0.004 
0.953 ±0.021 

1.066 = 0.032 


— — 1.003 ±0.011 

— — 0.972 ± 0.004 

— — 0.961 ±0.023 

l).97« ± 0.040 


Xtotal/DOF 
Xfkt/Nbin 
XSciFi/Nbin 
XSciBar/Nbin 


687.2 / 585 

85.4 / 80 
335.6 / 286 
266.1 / 239 


46.8 / 73 328.7 / 273 253.3 / 228 
47.7 / 80 — — 

— 328.7 / 286 - 

— — 253.3 / 239 



TABLE XV: The error matrix 
root of error matrix (sign [Mij] ■ 
unit of %. 



for /i and iinQE- The square 
y^pWiJf) is shown here in the 





fl 


/2 


fs 


/5 


fe 


fv 


/8 


RnQE 


/l 


43.86 


-3.16 


7.28 


-2.21 


-0.76 


-3.48 


0.81 


-8.62 




-3.16 


7.51 


1.97 


1.90 


0.62 


1.29 


2.43 


-5.68 


/3 


7.28 


1.97 


6.00 


3.38 


1.63 


3.44 


1.71 


-2.99 


h 


-2.21 


1.90 


3.38 


4.04 


-1.86 


4.53 


2.20 


1.65 


fe 


-0.76 


0.62 


1.63 


-1.86 


5.28 


-5.85 


5.11 


0.94 


fr 


-3.48 


1.29 


3.44 


4.53 


-5.85 


13.67 


-10.14 


4.09 


fs 


0.81 


2.43 


1.71 


2.20 


5.11 


-10.14 


18.35 


-11.77 


RnQE 


-8.62 


-5.68 


-2.99 


1.65 


0.94 


4.09 


-11.77 


20.30 




TABLE XVL SK event reduction summary. 



Reduction step 


K2K-I K2K-II 


|AT| < 500Aiscc, 


107892 


470469 


no pre- activity 






total number of p.e. within 


36560 


16623 


300 n sec timing window 






>200(K2K-I),94(K2K-II) 






Fully contained event 


153 


99 


flasher cuts 


97 


88 


visible Energy >30MeV 


95 


85 


fiducial volume cut 


56 


59 


|AT| = -0.2 - 1.3/xsec 


55 


57 




-500-400-300-200-100 100 200 300 400 500 
AT us 

FIG. 34: The AT distribution at each reduction step. Clear, 
hatched and shaded histograms are after pre-activity cut, to- 
tal p.e. cut, and fiducial volume cut, respectively. 
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TABLE XVII: SK event summary. For oscillated expecta- 
tions, sin^26l = 1 and Am^ = 2.8 x lO^^eV^ are assumed. 





K2K-I 


K2K-II 




data 


expected 


data 


expected 






w/o osc. 


w/ osc. 




w/o osc. 


w/ osc. 


Fully contained 


55 


80.8 


54.8 


57 


77.3 


52.4 


1-ring 


33 


51.0 


31.1 


34 


49.7 


30.5 


/i-like 


30 


47.1 


27.7 


28 


45.2 


26.7 


e-like 


3 


3.9 


3.4 


6 


4.5 


3.8 


multi-ring 


22 


29.8 


23.7 


23 


27.6 


21.9 



in 
c 
in 

« 

c 

0) 

>15 

UJ 



10 



n25ns 



5 




500 



1000 

AT (ns) 



ing structure of the beam can be clearly seen in the AT 
distribution if finer bins are used as in Fig. [35l Fig. 
shows the event rate as a function of POT. A KS-test 
performed against the assumption that the rate is pro- 
portional to POT gives a probability of 79%. 




10 20 30 40 50 60 70 80 90 100 

POT vs events CT *10" POT 



FIG. 36: Event rate as a function of POT. The KS-test prob- 
ability to observe our event rate under the assumption that 
the event rate is proportional to POT is 79%. 



FIG. 35: The AT distribution for fully contained events. The 
nine micro-bunch structure present in the beam is clearly seen. 



6. The events are selected to come from the 22.5 kilo- 
ton fiducial volume by requiring the reconstructed 
vertex position be at least 2m away from the inner 
detector wall. 

Tab . IXVTl shows the reduction summary for K2K-I and 
K2K-IL The efficiency for these cuts are 77.2% for K2K-I 
and 77.9% for K2K-II. The majority of the inefficiency is 
due to NC interactions which are not selected by these 
criteria. In total, 112 accelerator produced, fully con- 
tained events, are observed in the SK fiducial volume, 
with 58 events reconstructed as 1-ring /j,-like. Tab. IXVTll 
summarizes the characterization of these events and the 
MC expectations with and without neutrino oscillation. 

Figure [34] shows the AT distribution at each reduction 
step. A clear peak at AT = is seen after the fidu- 
cial volume cut. Three off-timing fully contained events 
are observed within ±500/zsec timing window which is 
consistent with the 2 expected background events from 
atmospheric neutrinos. In addition, the nine-bunch tim- 



The energy distribution of the events is compared 
against expectation in several ways. Fig. [37] shows the 
visible energy distribution, which is estimated from the 
energy deposit in the inner detector for all of the fully 
contained fiducial volume events. In this figure, the ob- 
served data is compared with the MC expectation based 
on the ND measurement without neutrino oscillation. 

Figure [38] shows the expected energy spectrum to- 
gether with the observation for the one-ring /i-like events. 
The expectation is normalized by the number of observed 
events (58). The neutrino energy is reconstructed us- 
ing the reconstructed muon momentum and the known 
beam-direction while assuming there was a QE interac- 
tion and ignoring the Fermi momentum. As can be seen, 
compared to the MC expectation there is a deficit of IR/i 
events in the low energy region, as is expected from the 
oscillation hypothesis. 



B. Systematic error from reconstruction at SK 

The systematic uncertainties for estimating NgK and 
the reconstructed neutrino energy in SK are evalu- 
ated using atmospheric neutrinos as a control sample. 
Tab. IXVIIII shows the systematic errors for Nsk- The 
dominant uncertainty for estimating Ngx comes from the 
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FIG. 37: The visible energy distribution for fully contained 
fiducial-volume events in SK. The closed circles are the ob- 
served data. The solid histogram is the MC expectation based 
on the ND measurement without neutrino oscillation, and the 
dashed one is the MC expectation with neutrino oscillation of 




e;^"^ GeV 

FIG. 38: The reconstructed Ei, distribution for the SK 1- 
ring /i-like sample. Points with error bars are data. The solid 
line is the expectation without oscillation. The histogram is 
normalized by the number of events observed (58). 



sin^26l = 1 and Am'' 



2.8 X IQ-^eV^ 



TABLE XVIII: Systematic errors for Nsk. 



vertex reconstruction. Since a cut is made on fiducial vol- 
ume a systematic shift in or out of this volume will either 
over or underestimate the number of events expected. It 
is evaluated comparing the number of events for atmo- 
spheric neutrino data with the MC expectation in the 
fiducial volume using two different vertex reconstruction 
programs. 



K2K-I K2K-II 


reduction 


<1% 


<1% 


fiducial volume cut 


2% 


2% 


decay electron background 


0.1% 


0.1% 


MC statistics 


0.6% 


0.6% 


Total 


3% 


3% 



Systematic errors due to the reconstruction algorithms 
themselves are also taken into account in the oscillation 
analysis. Systematic errors due to reconstruction are 
shown in table IXIXI Uncertainties coming from the ring 
counting and particle identification are evaluated by com- 
paring the likelihood distributions for data and MC, and 
varying the selection criteria. Figure [551 and HDl show the 
ring counting and particle identification likelihood dis- 
tributions of atmospheric neutrino data compared with 
the MC expectations in SK-II. The MC expectations re- 
produce the data well. The uncertainty for the energy 
scale are also estimated by using cosmic ray muons, the 
TT^ invariant mass and decay electrons. The energy scale 
uncertainty at SK is estimated to be 2.0% for K2K-I and 
2.1% for K2K-II. 



TABLE XIX: Systematic errors for reconstructed neutrino 
energy spectrum. The errors are shown in %, and the five 
columns refer to different bins of neutrino energy, as shown 
in the table in units of GeV. 



K2K-I (GeV) 


0-0.5 


0.5-1.0 


1.0-1.5 


1.5-2.0 


2.0- 


ring counting 


3.4% 


2.7% 


3.0% 


4.5% 


4.5% 


particle ID 


0.9% 


0.3% 


0.5% 


0.4% 


0.4% 


vertex 


2.0% 


2.0% 


2.0% 


2.0% 


2.0% 


total 


4.1% 


3.4% 


3.6% 


4.9% 


4.9% 




K2K-II (GeV) 


0-0.5 


0.5-1.0 


1.0-1.5 


1.5-2.0 


2.0- 


ring counting 


5.3% 


4.1% 


3.7% 


3.8% 


3.8% 


particle ID 


2.6% 


0.4% 


0.3% 


0.6% 


0.6% 


vertex 


2.0% 


2.0% 


2.0% 


2.0% 


2.0% 


total 


6.2% 


4.6% 


4.2% 


4.3% 


4.3% 
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IX. NEUTRINO OSCILLATION ANALYSIS 

A. Oscillation analysis method 

A two-flavor neutrino oscillation analysis is performed 
based on a comparison between the observation and the 
expectation by use of a maximum-likelihood method. 
The signatures of neutrinos oscillating from i/^ to Vr 
are both a reduction in the total number of observed 
neutrino events and a distortion in the neutrino energy 
spectrum. Thus, the likelihood function is defined as the 
product of the likelihoods for the observed number of 
events in the SK fiducial volume (/3noim) and the shape 
of the E'^j^'^ spectrum (£shapo)- In addition, the system- 
atic uncertainties are also treated as fitting parameters 
in the likelihood. They are included in a constraint like- 
lihood term (>Csyst) where they are varied within their 
uncertainties, thus modifying the expectation. The total 
likelihood function is then defined as: 



FIG. 39: Ring counting likelihood distribution for SK-II at- 
mospheric fully contained neutrino events. Closed circles are 
data and the histogram is MC expectation normalized by live 
time assuming neutrino oscillation at atmospheric best-fit pa- 
rameters. Events with Likelihood < are assigned to be one- 
ring. 
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FIG. 40: Particle identification likelihood distribution for SK- 
II atmospheric fully contained one-ring events. Closed circles 
are data and the histogram is MC expectation normalized by 
live time assuming neutrino oscillation at atmospheric best-fit 
parameters. Events with Likelihood < are assigned to be 
e-like. 



X C. 



shape 



X C 



syst- 



(21) 



The oscillation parameters, Am^ and sin^ 26', are ob- 
tained by maximizing the likelihood function. One- 
hundred twelve FC events are used in £norm and fifty 
eight FC IR/J, events are used for £shapc, respectively. 
The systematic parameters in the likelihood consist of 
the neutrino energy spectrum at the near detector site, 
the F/N flux ratio, the neutrino- nucleus cross-section, 
the efficiency and the energy scale of SK, and the overall 
normalization. 



Prediction of the number of events and the 
energy spectrum at Super-Kamiokande 



Number of neutrino events 



The expected number of neutrino events in SK {N^^) 
is derived by extrapolating the measured number of inter- 
actions in the IKT {N^^^ , calculated in Eq.il) with the 
ratio of the expected neutrino event rate per unit mass, 
pSK I piWT ^ Taking into account the difference of fiducial 
mass {M) and the number of protons on target (POT) 
used in the analysis for IKT and SK, is written as: 



IKT 



M' 



SK 



POT 



SK 



3IKT MiKT poT^^'^ 



■a., (22) 



where superscripts "IKT" and "SK" denote the variable 
for SK and IKT, respectively, and Ci,^ is the correction 
factor for the difference of the electron neutrino contam- 
ination in the neutrino beam at IKT and SK. The value 
of Cj/ is estimated to be 0.996 with the MC simulation. 
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The expected event rate at each detector, p, is calcu- 
lated from the neutrino flux <&, the neutrino- water in- 
teraction cross-section cr, and the detector efficiency e 
estimated with the MC simulation: 



p = J dE,^E,)-a{E,)-e{E,). 

In order to account for the systematic uncertainty, we 
classify the neutrino interaction into three categories: 
CC-QE, CC-nonQE and NC. The event rate is calcu- 
lated separately for each of the three interaction types 
and then summed. The neutrino flux at SK, <i>^^, is es- 
timated from the F/N flux ratio i?^/^and the measured 
ND spectrum $^°: 



We estimate the expected number of events without 
neutrino oscillation while incorporating all of the known 
systematic uncertainties by use of a MC technique. Many 
sets of the systematic parameters are randomly generated 
with proper consideration of their correlation. For each 
systematic parameter set, is calculated using Equa- 
tion (22), setting oscillation parameters to zero. The 
number of EC events without neutrino oscillation is es- 
timated to be 158.ll;g;g. This technique allows us to 
determine the contributions from the individual system- 
atics to the the total error by selectively including only 
some errors during the generation. We find that the dom- 
inant error sources are the fiducial volume uncertainties 
in IKT and SK (+4]^°) and the F/N ratio (+2:9!) ■ 



^SK ^ i?F/N(^^) . $ND(£;^) . _ P{E,;Am^,sm^ 29)), 

where P{E^; Am^ ,sm^ 29) is the neutrino oscillation 
probability given by Eq. ([1]). 

The uncertainties of and their correlation are ob- 
tained from the ND analysis as shown in Tab. lXV] Those 
for R^/^ are derived from the HARP 7r+ measurement 
and the beam MC simulation, and are summarized in 
Tab. lYl In order to be insensitive to the absolute cross- 
section uncertainty, we incorporate the uncertainties in 
neutrino-nucleus cross-section as the cross-section ratio 
relative to CC-QE interactions. The uncertainty of the 
CC-nonQE/CC-QE cross-section ratio is taken from the 
ND measurements. For the NC/CC-QE cross-section 
ratio, we assign 11% uncertainty to NC single 7r° pro- 
duction based on the measurement with the IKT [67j . 
The other NC interaction modes are assigned 30% un- 
certainty based on past experiments (68j . Taking into 
account the detection efficiency in SK, 15% is assigned 
as the net uncertainty on NC/CC-QE ratio. The uncer- 
tainties from event reconstruction at SK are summarized 
in Tab. IXVIIII The uncertainty of the overall normaliza- 
tion of the number of events in each period is estimated 
from the fiducial mass error of IKT and SK and the un- 
certainty in the difference of the number of protons on 
target used for the analysis. In total, the normalization 
error is estimated to be ±5.1% for both lb and II periods. 

For the period la, the beam configuration is different 
from other periods and we do not have a ND measure- 
ment of the energy spectrum. We employ different treat- 
ment of the systematic errors for this period. All the 
uncertainties, including those from the F/N flux ratio, 
energy spectrum and cross-section, are incorporated into 
the error of the single normalization parameter. The 
total normalization uncertainty for the la period is es- 
timated to be -|-9.0/-9.8%. The largest contribution to 
this uncertainty is from the energy spectrum, which is es- 
timated with the HARP vr"*" measurement and the beam 
MC simulation and gives -|-5.8%/-7.0%. Other contri- 
butions come from the F/N flux ratio (±4.3%) and the 
fiducial volume uncertainties of IKT (±4.3%) and SK 
(±3.0%). 



2. Reconstructed neutrino energy spectrum 

The expected spectrum shape of the reconstructed 
neutrino energy at SK, (f>^^p{El'^^) , is calculated as: 

</>SK ^ JdE,- $ . ^(^E,) . ef^^iE,) ■ T{E,-Ein, 

(23) 

where e^^^ is the detection efficiency for IR/i events in 
SK and r{Ei,; E'J"^) is the probability of reconstructing 
an event with true energy E^, as Both of them 

are estimated with the MC simulation. In the actual 
procedure, the Ei, and E^'^'^ are binned with an interval 
of 50 MeV, and hence the integral over the true neutrino 
energy is replaced by a summation over true energy bins. 

The uncertainties from the neutrino energy spectrum 
at the ND, the F/N fiux ratio, and the cross-section ra- 
tios are incorporated as described above. The uncertain- 
ties from IR^ event reconstruction at SK are shown in 
Tab. IXIXI The energy scale uncertainty in SK is 2.1% 
for SK-I and 2.0% for SK-II, respectively, as described in 
section 

The expected i?"'^ spectrum shape for null oscillation 
case and its error are estimated using the same technique 
as the number of events and shown in Fig. 311 The height 
of the box represents the size of estimated error in each 
bin. The contribution of each systematic uncertainty is 
estimated by turning each uncertainty on exclusively one 
by one, as shown in Fig. [151 We find that the error on 
the spectrum shape is dominated by the SK energy scale. 



C. Definition of likeliliood 

1. Normalization term 

The normalization term, £norm, is defined as the Pois- 
son probability to observe A^obs events when the expected 
number of events is N^xv'- 
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FIG. 41: Expected reconstructed neutrino energy spectrum 
shape in the case of null oscillation. Height of boxes indicate 
the size of error. 



^^fol 



-N„: 



(24) 



In order to account for the difference of the experimen- 
tal configuration, the expectation for each experimental 
period is separately calculated using Eq. ((22|) and then 
summed as: 




FIG. 42: Contribution of each systematic error to the recon- 
TVox — N^^ + N^^ + N^^ (25) structed neutrino energy spectrum. Vertical axis is relative 

cxp exp exp error of the spectrum. Source of uncertainty is indicated in 

each plot. Blank and filled bars represent the sizes of the total 
error and the contribution from the source being considered, 
2. Energy spectrum shape term respectively. 



The energy spectrum shape term is defined as the prod- 
uct of the probability for each IR/i event to be observed 
at reconstructed neutrino energy E'J"^. We use the ex- 
pected neutrino energy spectrum, given in Eq. (|23p. as 
the probability density function. The probability den- 
sity function is separately defined for each experimental 
period: 

/:.hapc = n 0cfp.ib(^^:?^;Am2,sin2 20) 

i=l 

X n ^t^p,u{El7;^rri',sm'2e) (26) 

where Nl^^ — 30 and Nl]^^ = 28 are the number of 
observed FC IR/i events for period lb and II, respectively. 
There is no IKfi event in the la run period. 



3. Systematic term 

The systematic parameters are treated as fitting pa- 
rameters, and are assumed to follow a Gaussian distribu- 
tion. They are constrained within their uncertainties by 
constraint terms expressed as: 



n exp(-A//(Af,r'A/,.), (27) 

where iVgyst is the number of parameter sets, A/^ repre- 
sents the deviations of the parameters from their nominal 
values and Mj is the error matrix for j-th set of param- 
eters. 
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D. Results 

The likelihood is maximized in the Am^ - sin^ 29 space 
and the best fit point within the physical region is found 
to be at (Am2,sin2 26i) = (2.8 x 10-3eV^1.0). The val- 
ues of all systematic parameters at the best fit point 
are within la of their estimated errors. At this point, 
the expected number of events is 107.2, which agrees 
well with the 112 observed within the statistical uncer- 
tainty. The observed El°^ distribution is shown in Fig. [43] 
together with both the expected distributions for the 
best-fit parameters, and the expectation without oscil- 
lations. The consistency between the observed and the 
best-fit El^"^ distributions is checked using a Kolmogorov- 
Smirnov (KS) test. For the best fit parameters, the KS 
probability is 37 %, while for the null oscillation hypothe- 
sis is 0.07 %. The observation agrees with the expectation 
of neutrino oscillation. The highest likelihood is found at 
(Am^, sin^ 20) = (2.6 x lO-^eV"^, 1.2), which is outside of 
the physical region. The probability that we would get 
sin^ 20 > 1.2 if the true parameters are at our best fit 
point is 26.2%, based on the virtual MC experiments. 
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TABLE XX: Summary of the null oscillation probability. 
Each row is classified by the likelihood term used, and each 
column represents the data set. 



K2K-I+II K2K-I only K2K-II only 

Shape + Norm. 0.0015% (4.3cr) 0.18% (S.lcr) 0.56% (2.8cr) 

Shape only 0.42% (2.9cr) TWo 5^2% 

Norm, only 0.06% (3.4o-) 0.6% 2.8% 



TABLE XXL Effect of each systematic uncertainty on the 
null oscillation probability. The numbers in the table are null 
oscillation probabilities when only the error written in the 
first column is turned on. 



Norm-only Shape-only Combined 



Stat, only 


0.01% 


0.22% 


0.0001% 


FD spectrum 


0.01% 


0.24% 


0.0002% 


uQE/QE, NC/CC 


0.01% 


0.23% 


0.0002% 


Far/Near 


0.02% 


0.23% 


0.0003% 






0.23% 


0.0002% 


Energy scale 




0.38% 


0.0002% 


Normalization 


0.03% 




0.0005% 


All errors 


0.06% 


0.42% 


0.0015% 



The null oscillation probability calculated separately 
for each sub-sample or each likelihood term is shown in 
Tab.KXl In addition. Tab. IXXll shows the effect of each 
systematic uncertainty on the null oscillation probability. 
The effect is tested by turning on the error source written 
in the first column in the table. As shown in the table, 
the dominant contributions to the probabilities for the 
normalization information are from the F/N fiux ratio 
and the normalization error, while the energy scale is 
the dominant error source for the probability with the 
El°^sha.pe information consistent with the results found 
using the MC test described in Sec. IIXB2I 

The allowed region of oscillation parameters are eval- 
uated based on the difference of log-likelihood between 
each point and the best fit point: 



FIG. 43: The reconstructed Ei, distribution for the 1-ring fi- 
like sample. Points with error bars are data. The solid line is 
the best fit spectrum with neutrino oscillation and the dashed 
line is the expectation without oscillation. These histograms 
are normalized by the number of events observed (58). 



The probability that the observations can be explained 
equally well by the no oscillation and by the oscillation 
hypotheses is estimated by computing the difference of 
log-likelihood between the null oscillation case and the 
best fit point with oscillation. The null oscillation prob- 
ability is calculated to be 0.0015 % (4.3o-). When only 
normalization (shape) information is used, the probabil- 
ity is 0.06% (0.42%). 



/ /-phys \ 

Aln£(Am2,sin2 26l) = In ti^^ 

= In £^^^ - In £( Am^ , sin^ 26) , 

(28) 

where Cf^^^ is the likelihood at the best-fit point and 
C{Am'^,sin^2e) is the likelihood at (Am^, sin^ 26*) with 
systematic parameters that maximize the likelihood at 
that point. 

The allowed regions in the neutrino oscillation param- 
eter space, corresponding to the 68%, 90% and 99% con- 
fidence levels (CL) are shown in Fig. They are de- 
fined as the contour lines with ln£ = ln£Pj!^^ — 1.37, 
—2.58 and —4.91, respectively. These regions are derived 
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TABLE XXII: Summary of the oscillation parameters at the 
best fit point for each fitting condition. 



All region 



Physical region 



All data shape + norm. 2.55x10"'^ 1.19 2.75x10 



shape only 2.77x10-^ L25 2.95x10 



1.00 



1.00 



K2K-I shape + norm. 2.77x10' 



K2K-II shape + norm. 2.36x10"'' 1.35 2.64x10 



1.08 2.89x10 



1.00 



1.00 



by using the two-dimensional Gaussian approximation 
from the maximum in the unphysical region (2l| . The 
90% C.L. contour crosses the sin^ 26 — 1 axis at Arn^ = 
1.9 and 3.5 X 10-3 eV^ Fi gure l45l shows the distributions 
of InCf^y^ - ln£(Am2,sin2 26) as a function of sin^ 29 
and Ato^, with a slice at either Am^ = 2.8 x lO^'^eV^ or 
sin^ 26 = 1.0. 

We also check the consistency of the fit results per- 
forming the analyses with only the normalization term 
or spectrum shape term, and with the K2K-I or K2K-II 
sub-samples separately. The fit results are summarized in 
Tab. lXXlTl There is no entry for the normalization term 
only, because the two parameters cannot be simultane- 
ously determined from only one number. The oscillation 
parameters allowed by the normalization and the spec- 
trum shape alone agree with each other, as shown in both 
Tab. IXXIII and Fig. [351 The allowed regions calculated 
with only K2K-I and K2K-II data are also consistent as 
shown in Tab. IXXlIl and Fig. gT] 

Finally, we compare our result with the parameters 
found by the measurement of atmospheric neutrino oscil- 
lation by the Super-Kamiokande collaboration 3- Fig- 
ure|48]shows the allowed regions of oscillation parameters 
found in this analysis together with the SK result. The 
K2K result is in good agreement with the parameters 
found using atmospheric neutrinos, thereby confirming 
the neutrino oscillation result reported by SK. 



X. SUMMARY 

Data taken by the K2K experiment between June 1999 
and November 2004 is used to observe and measure the 
parameters of neutrino oscillation using an accelerator- 
produced neutrino beam. The K2K experiment is the 
first long-baseline neutrino experiment to operate at a 
distance scale of hundreds of kilometers. The neutri- 
nos are measured first by near detectors located approx- 
imately 300 meters from the proton target, and then by 
the Super-Kamiokande detector 250 km away. The near 
detector complex consists of a 1 kiloton water Cherenkov 
detector, and a fine grained detector system. The energy 
spectrum and flux normalization measured at the near 
detectors are used to predict the signal seen at Super-K. 
The results found are consistent with the neutrino os- 
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FIG. 44: Allowed regions of oscillation parameters. Three 
contours correspond to the 68% (dotted line), 90% (solid line) 
and 99% (dashed line) CL. allowed regions, respectively. 
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sin 29 



FIG. 45: InLgJS'x" - ln/:(Am^sin^ 29) distribution as a func- 
tion of Am^ (left) and sin^ 29 (right), sin^ 29 is set to be 1.0 
in the left-hand figure and Am^ is set to be 2.8 x IQ-^eV^ 
in the right-hand figure. Three horizontal lines correspond 
to the 68%, 90% and 99% CL interval from the bottom one, 
respectively. 



dilation parameters previously measured by the Super- 
Kamiokande collaboration using atmospheric neutrinos. 

One hundred and twelve beam-originated neutrino 
events are observed in the fiducial volume of Super- 
Kamiokande with an expectation of 158.llg5 events 
without oscillation. The spectrum distortion expected 
from oscillation is also seen in fifty-eight single-ring 
muon-like events which have had their energy recon- 
structed. A likelihood analysis was performed and 
the probability that the observations are explained by 
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Number of events only 



Spectrum shape only 




0.25 0.5 0.75 1 
Sin^(2e) 



FIG. 46: Allowed region of oscillation parameters evaluated 
with the number of events only (left) and the El°'^ spectrum 
shape only (right). Both information allow the consistent 
region on the parameters space. 



K2K-I only 



K2K-II only 




0.25 0.5 0.75 1 
Sin^(2e) 



0.25 0.5 0.75 1 
Sin^(2e) 



FIG. 47: Allowed region of oscillation parameters evaluated 
with partial data of K2K-I-only (left)/K2K-II-only (right). 
Both data allow the consistent region on the parameter space. 
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a statistical fluctuation with no neutrino oscillation is 
0.0015% (4.3(t). In a two flavor oscillation scenario, the 
allowed Am^ region at sin^ 26* = 1 is between 1.9 and 
3.5 X 10-^ eV^ at the 90 % C.L. with a best-fit value of 
2.8 X 10-^ eV2. 



FIG. 48: Comparison of K2K results with the SK atmo- 
spheric neutrino measurement Dotted, solid, dashed and 
dash-dotted lines represent 68%, 90%, 99% C.L. allowed re- 
gions of K2K and 90% C.L. allowed region from SK atmo- 
spheric neutrino, respectively. 
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